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r~| ' A practical approach is presented which allows the use of a non-invariant regularization 

P^, scheme for the computation of quantum corrections in perturbative quantum field theory. 

CIhI The theoretical control of algebraic renormalization over non-invariant counterterms is 

^ I translated into a practical computational method. We provide a detailed introduction 

into the handling of the Slavnov- Taylor and Ward-Takahashi identities in the Standard 
Model both in the conventional and the background gauge. Explicit examples for their 
practical derivation are presented. After a brief introduction into the Quantum Action 
5^ I Principle the conventional algebraic method which allows for the restoration of the func- 

tional identities is discussed. The main point of our approach is the optimization of this 
procedure which results in an enormous reduction of the calculational effort. The coun- 
terterms which have to be computed are universal in the sense that they are independent 
of the regularization scheme. The method is explicitly illustrated for two processes of 
phenomenological interest: QCD corrections to the decay of the Higgs boson into two 
photons and two-loop electroweak corrections to the process B — > ^^7. 
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1 Introduction 



A regularization method which respects all symmetries of the Standard Model (SM) [|l| 
does not exist. The popular and powerful method of Dimensional Regularization |^ is at 
least an invariant scheme for QCD. In the electroweak sector, however, the coupling to chiral 
fermions introduces the well-known 75 problem. One also has to face additional technicalities 
due to evanescent operators in the effective field theory approach. It is well-known that 
in the framework of Dimensional Regularization only the 't Hooft-Veltman-Breitenlohner- 
Maison scheme Q for 75 is shown to be consistent to all orders. The so-called naive 
dimensional scheme (with an anticommuting 75) does not reproduce the chiral anomaly and 
is not consistent to all orders. For specific examples it leads to correct results at the lowest 
orders in perturbation theory. Nevertheless, it seems desirable to have a powerful practical 
alternative even in the SM, at least for cross-checks, as suggested controversies in the past 
suggest (see, e.g., 0). Moreover, the impressive experimental precision mainly reached at the 
electron positron colliders LEP and SLC and the proton anti-proton collider TEVATRON has 
made it mandatory to evaluate specific two- or even three-loop contributions to observables 
where the inconsistencies of the well-known "naive dimensional" scheme are unavoidable. 

Going beyond the SM, it is well-known that Dimensional Regularization breaks the Ward 
identities of supersymmetry. However, one very often prefers to keep the dimensional scheme 
for the practical calculations, also beyond the SM, in order to take advantage of already well- 
developed computer tools 0. Thus, one needs a practical procedure to restore the Ward 
identities of supersymmetry in the final step of the renormalization procedure. 

From the principal point of view, the calculation of higher-loop contributions in pertur- 
bative quantum field theories is a well-understood issue. The axioms of relativistic quantum 
field theory, such as causality and Poincare invariance, fix the matrix elements completely 
to all orders up to a limited number of free constants. They have to be determined by 
renormalization conditions. These free constants correspond to a renormalization ambiguity 
for coinciding points in the definition of time-ordered products of operator-valued distribu- 
tions The main question is whether the renormalization ambiguity can be fixed in such 
a way that the time-ordered products fulfill the symmetry constraints. The question behind 
this is the compatibility of the symmetries of the classical Lagrangian with quantization. 

Here the method of algebraic renormalization offers a complete theoretical answer: In 
general, the subtraction of ultra-violet divergences in quantum field theories leads to non- 
invariant Green functions, which means that the regularization scheme and the subsequent 
renormalization do not respect the symmetries of the theory like supersymmetry or local 
gauge symmetries. As we mentioned above. Dimensional Regularization preserves gauge 
symmetries (up to the 75 problem) but breaks supersymmetry. 

The Quantum Action Principle [^] tells us that the breaking terms are local at the lowest 
non-vanishing order. This fact provides a possible path for the construction of invariant 
Green functions, independent of the regularization scheme. One introduces, order by order, 
finite non-invariant local counterterms which restore the symmetry relations (provided there 
are no anomalies) [Q. Thus, one can in principle show that in anomaly- free theories the 
local renormalization ambiguity (which is not fixed by the axioms of relativistic quantum 
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field theory) can always be used in such a way that the perturbative S-matrix enjoys all 
symmetry properties of the classical theory (for a review, see |1T0[| ). 

Although the method of algebraic renormalization is intensively used as a tool for proving 



renormalizability of various models [|I0|, its full value has not yet been widely appreciated by 
the practitioners. Indeed, the theoretical understanding of algebraic renormalization does 
not lead automatically to a practical advice for higher-loop calculations. One could even 
expect that such an algebraic renormalization scheme becomes very complicated at higher 
orders. It is one of the main purposes of this paper to provide theoretical procedures which 
minimize the additional efforts for the restoration of the symmetries and to demonstrate the 
efficiency of the combined method in some examples of phenomenological interest. 

However, two obvious practical complications of algebraic renormalization have to be 
taken into account: 



(a) The constraints introduced by the symmetry connect various Green functions. Thus, 
for the construction of the non-invariant counterterms corresponding to a specific Green 
function one also has to compute the various other Green functions involved in the 
identities. 

(b) In the computation of higher-loop contributions one also has to analyze identities from 
lower orders which constrain the non-invariant counterterms. 



These disadvantages can be significantly reduced: 

(1) First, one should state that many identities are not relevant if one is interested in one 
specific Green function only and if the corresponding breaking terms can be compen- 
sated by the other Green functions in the given identity alone. 

(2) In the case of local gauge symmetries, the structure of the relevant identities can be 



considerably simplified by using the background field gauge |TT[]. In a conventional 
gauge there is a large number of non-linear Slavnov- Taylor identities. In the back- 
ground field gauge some of them get replaced by linear Ward-Takahashi identities like 
in QED. 



(3) We have some well-known theoretical constraints [|TU]: the Quantum Action Principle 
tells us that the breaking terms are local at the lowest non-vanishing order and thus 
are removable by counterterms if there is no anomaly. Furthermore, the algebraic 
consistency conditions heavily constrain the structure of the breaking terms. 

(4) Finally, the most important simplification we want to present in this paper is the 
following: the number of breaking terms one has to calculate in addition can essentially 
be reduced to the ones which correspond to finite Green functions. This can be achieved 
by using a specific zero-momentum subtraction procedure. 

In this paper we want to discuss these different ingredients from a practical point of view 
and offer an algorithmic strategy for practical algebraic renormalization. As illustrating 
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examples for our combined algebraic method we have chosen two processes of phenomeno- 
logical interest, namely the two-loop contributions to B ^ X^'y and to H ^ 77. The 
important extensions of these techniques to supersymmetric examples will be presented in a 
forthcoming paper. 

As mentioned above, the proposed procedure based on algebraic renormalization is not 
restricted to a specific class of regularization schemes. Once the structure of the local 
breaking terms are under control, one can choose the most practical regularization scheme 
for the specific case under consideration. 

In the following we also use the method of Analytic Regularization in one of our illustrat- 
ing examples. This choice is guided by the fact that this scheme enjoys the property of mass 
independence like the minimal subtraction (MS) @, ^ or the modified minimal subtraction 
(MS) scheme of the Dimensional Regularization. 

The delicate infra-red problem is another important task. As mentioned above, the 
method includes zero-momentum subtractions which heavily rely on the regularity properties 



of the Green functions at zero momentum [|T^]. Here we mention the necessary modifications 
in massless theories. 

The paper is organized as follows: 

In Section |] we recall the fundamental symmetry constraints of the SM namely the 
Slavnov- Taylor and the Ward-Takahashi identities. The main idea of this chapter is to 
collect all technical ingredients which are necessary to derive the symmetry constraints for 
a specific process in the SM. 

In the first part of Section ^ we discuss the practical consequences of two further ingredi- 
ents of the algebraic renormalization, namely the Quantum Action Principle and the Wess- 
Zumino consistency conditions, particularly within the background field method (BFM). 
Then we propose our main procedure to remove the breaking terms in the specific symmetry 
identities. The various practical steps are presented in an algorithmic form. 

In Section § we illustrate our practical algebraic renormalization scheme in the two-loop 
calculation of the decay H —>■ 77, which is one of the promising channels for the discovery 
of the Higgs boson with a mass of around 120 GeV. 

In Section |^ the analysis of the electroweak corrections to the decay 6 — > S7 is presented. 

In the Appendices some auxiliary technical and theoretical information used in Sections ^ 
and 1^ are offered to the reader. In particular those parts of the SM Lagrangian in the 
background field gauge which are absent in the literature are given in Appendix ^ In 
Appendix an explicit example on how in practice the Slavnov- Taylor identities are derived 
is discussed. In Appendix y we analyze the triangular structure of the counterterms further. 
This analysis allows to restore the identities in a step-by-step procedure. 



2 Slavnov- Taylor and Ward-Takahashi identities 

The main tools for algebraic renormalization are the Slavnov- Taylor (STI) and Ward- 
Takahashi identities (WTI). In this Section it is shown how the complete set of identities 
corresponding to a specific process are derived from their general form and how it is possible 
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to disentangle the contributions coming from QCD and electroweak radiative corrections. 

At this point a word concerning the notation is in order. A generic field is denoted 
by (p. $ stands for scalar matter fields, i.e. Goldstone (G^, G°) and Higgs bosons (H). 
Fermionic fields, respectively their conjugates are represented by and ip. A generic gauge 
boson field is denoted by V/^ and the ghost and the anti-ghost fields by c and c, respectively. 
The symbols and c** are used to denote gluon fields and the corresponding ghosts in the 
adjoint representation of the Lie algebra su{3). The background fields are marked with a hat 
in order to distinguish them from their quantum counterparts. Qi, respectively Qg-, denotes 
the electric charge of a quark qi. 

Let us also introduce three different types of effective actions which will be used in the 
following. The Green functions F are regularized and renormalized. The Green functions F 
are subtracted using Taylor expansion (see Section p.3| ) . Finally, F denotes the renormalized 
symmetric Green functions, which satisfy the relevant WTIs and STIs. 

A complete explanation of the conventions, quantum numbers and symmetry transfor- 
mations of the fields is provided in Appendix ^ 



2.1 Conventional gauge fixing 

In this Section the general form of the STI in the conventional 't Hooft gauge fixing is 



presented. Thereby we follow the so-called Zinn- Justin formalism ||I4 . 

Let us consider the Gell-Man-Low formula for one-particle irreducible (truncated) Green 
functions (IPI) 

= (T(0°(xl)...0:(x„))e-/^'-^-(-))l^^ (2.1) 

where the superscript "o" recalls the free fields. The Fourier transformed Green functions are 
denoted by T^^ ^^{pi, . . . ,p„) where pi, . . . ,pn are the incoming momenta^. The definition 
of T^^,,,^^{pi, . . . ,pn) in terms of time-ordered products of free fields, (pl - . . 0°, and vertices 
of the interacting Lagrangian, Cint, requires a regularization and a subtraction prescription. 
In this Section we do not rely on a specific scheme, but only on general features of the 
renormalization theory such as the Quantum Action Principle (QAP) (see Section |3]l|) and 
the Zimmermann identities |15|. 



To handle the complete set of Green functions, it is very useful to collect them into a 
generating functional 



ITM = E / n Y.P^ Mpi) ■ ■ ■ MPn)r^,...^„{pu ...,Pn). (2.2) 

n=0 \j=0 / \ k / 

In perturbation theory r<^-^,..0^(pi, . . . ,pn) is a formal power series in h. In the following, we 
will adopt the notation T^^\^^ to indicate the m-loop contribution to the Green function 

^Here and in the following momentum conservation is assumed, i.e. Pi = 0- 
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(a) (b) 

Figure 1: All momenta are considered as incoming. In the Green functions they are 

assigned to the corresponding fields starting from the right. The momentum of the most left field 
is determined via momentum conservation. Exemplary T^^^^{p) and ^^^^.^^^(g,^) are pictured 
in (a) and (b), respectively. 



^(j)i...(l>n- terms of F each single Green function of the form ( |2.1| ) is obtained by means of 
functional derivatives 



501 (Pl) . ..6(j)n{Pr. 



= r<^l...<An(Pl • • -Pn), (2.3) 

<f)=0 



where (f){p) denotes the Fourier transform of 0(x). In Fig. |l| our conventions concerning the 
external momenta can be found. The Green functions of Eq. ( |2.1| ) exhaust all the possible 
amplitudes involved in the S-matrix computation, but they do not cover the complete set 
of Green functions needed for the renormalization of the theory. Indeed, due to the non- 
linearity of the Becchi-Rouet-Stora-Tyutin (BRST) transformations 0, the renormalization 
of some composite operators (namely s0j where (pi is a generic field of the SM and s is the 
BRST generator) is necessary. This is usually done by adding the composite operators scpi 
coupled to BRST invariant external sources 0* to the classical action 



where Cmv is the gauge invariant Lagrangian of the SM (see |1T6| , \V7[ p!q| and remarks in 
Appendix ^) and studying the renormalization of C. For our purposes we only introduce the 
BRST sources (also called anti-fields) for non-linear transformations as proposed by Zinn- 
Justin As a remark we mention that in the Batalin-Vilkovisky anti- field formalism |^ 



the BRST sources are also introduced for linear BRST transformations. The advantage is 
that all the gauge fields occur on the same footing. However, they are neither necessary for 
our practical purposes nor for proving the renormalization of the SM. 

The quantization of the theory can only be achieved by introducing a suitable gauge 
fixing Cgf and the corresponding Faddeev-Popov terms £$n 

To = / d^x (ciNv + E + ^GF + Un] . (2.5) 
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Both Cgf and £$n break the local gauge invariance leaving the theory invariant under 
the BRST transformations. The ERST symmetry is crucial for proving the unitarity 
of the S-matrix and the gauge independence of physical observables. Therefore it must 
be implemented to all orders. For this purpose we establish the corresponding STI in the 
functional form 
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(2.6) 



where the notation A'^B^ = A~^B~ + A~B^ has been used, sw and cw denote the sine and 
cosine of the Weinberg angle 9w and are the so-called the Nakanishi-Lautrup multiplier^ 
The sum in the last line of Eq. ( |2.6| ) includes the left-handed doublets and the right-handed 
singlets. For the BRST source fields no Weinberg-rotation has been introduced. We stress 
that this formula represents the complete nonlinear STI to all orders. The first two and the 
last term correspond to the linear BRST variation of the U{1) abelian gauge field and the 
BRST transformations of the anti-ghost fields. Note that the STI of the form ( ^.6| ) contains 
the complete information of the BRST symmetry and the equation of motion p|, |14[. 

In the form of Eq. ( p.6|) the STIs are independent from the gauge fixingf] In order to 
specify the gauge fixing, we introduce the equation of motion for the h fields corresponding 
to the various gauge fields in the SM 



5T_ 

6ba 



J-,(V,<l>)+ea&o 



(2.7) 



where JF„ (a = A, Z, W, g) are the gauge fixing functions. (« = A, Z, g) are the 
corresponding gauge parameters. In the case of the background gauge fixing the functions 
JFq, are explicitly given in the formula (A. 2) of the Appendix. 

Considering a specific process, one first has to single out the complete set of relevant 
identities by using a functional derivative (as in Eq. (|2.3|)). With relevant set we mean the 



^In practical calculations they can be eliminated (in the case of linear gauge fixing) by a Gaussian 
integration. 



•^Note that we do not have to modify Eq. (2.6) if the gauge fixing is changed from the conventional 't Hooft 
gauge (see Section 2.2) to the background gauge which is used in Section 2.4. However, in order to control 



the dependence of the Green functions on the background fields some new terms are conventionally added 
to the STIs. They implement the equation of motion for the background fields and they are studied in the 



Section 2.3 
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set of identities which is closed under renormahzationQ. This means that the finite parts 
of a Green function appearing in a given identity is fixed by other identities of the set or 
by renormahzation conditions. In practical calculations usually not all identities are really 
necessary since they might be automatically preserved by reasonable regularization schemes 
at lower orders. An identity can also decouple from the others, because it only contains 
Green functions which do not infiuence the breaking terms of the other identities. The latter 
point will be discussed in Section ||. 

The most convenient procedure to deduce the complete set is the following: (i) consider 
the amplitudes involved in the physical process; (m) derive the identities for those amplitudes; 
[iii) from each identity single out the new (superficial divergent) Green functions which are 
not involved in the physical process; (iv) derive the identities for these new Green functions. 
The procedure stops when the new identities involve only new finite Green functions and 
no other divergent quantities. Finally we have to underline that supplementary constraints 
such as the Faddeev-Popov equations can lead to relevant identities on Green functions which 
avoid the use of a new STL E.g., in the case of the two-point functions no derivative w.r.t. 
ba has to be considered. 

In order to obtain a meaningful expression the following two simple rules have to taken 
into account: 

1. Green functions with a positive or negative Faddeev-Popov ghost charge vanish as it is 
conserved. Thus, in order to extract non-zero identities, it is necessary to differentiate 
the expression 5(r) = 0, which carries ghost charge w.r.t. one ghost field also 
having ghost charge +1. It is also possible to differentiate w.r.t. two ghost fields and 
one anti- field (carrying ghost charge —1). The only exception to this rule is the case of 
anti-fields for the ghosts. They carry two Faddeev-Popov ghost charges and, therefore, 
these charges must be compensated with three ghost fields. 



2. Identities for the Green functions are obtained by taking derivatives of the STI ( |2.6| ) 
w.r.t. fields and external sources. Clearly, they are non-vanishing only if Lorentz 
invariance is respected. 

The derivation of the complete set of non-trivial identities is guided by the following rules: 

3. If we are interested in identities involving several gauge bosons one has to differentiate 
5(r) = w.r.t. the set of fields where one of the gauge bosons is replaced by the 
corresponding ghost field Cj. The reason for this is that the linear part of the BRST 
transformations of a gauge field is proportional to the corresponding ghost: sV^ = 
di,c+ .. .. 

4. For Green functions which contain ghost fields a new rule is needed. One ghost field 
must be replaced by two ghost fields. In fact, the BRST transformation of the ghost 
fields is non-linear scj = ^fijkC^c'' where fijk are the structure constants of the gauge 
group. In the case of ghost two-point functions this is not necessary because we do not 
acquire any new constraints on them from this rule (see also Appendix 0). 



*Up to additional relations which get eventually introduced by normalization conditions. 
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5. In the identities derived with the help of rules ^ and ^ different Green functions occur. 
The ones which still involve gauge bosons or ghosts are constrained further by identities 
which one may derive as described above. 

Since the STIs will be our main tools in the context of algebraic renormalization, we 
want to consider their derivation from Eq. (|2.6| ) in more detail. Recall that both r[0] and 
iS(r)[0] are integrated functional of the fields 0. Thus it is possible to apply the rules of 
functional derivatives (see, e.g., Ref. Section 6-2-2). Taking the functional derivatives of 



r and setting afterwards all fields to zero generates a single Green function r^j,..,^^. On the 
other hand, the functional derivatives of 5(r)[0] generate a single STI (again after setting 
the fields to zero after differentiation). Note that in the expression for iS(r)[0] already some 
functional derivatives are present which must be interpreted as functionals of the form ^[0]- 
The use of rules for taking the derivative of products enables us to distribute the functional 
derivatives to the individual expressions in S (F) [0] and to set all fields to zero afterwards. 

From the technical point of view the only detail to be clarified is the dependence on the 
space-time coordinate, respectively, the momenta of each single STI. The presence of the 
integral over the space-time in Eq. ( |2.6| ) and the conservation of the momentum flow of the 
Green functions guarantees that no momentum integration is left. Thus the STI can be 
expressed as a sum of products of Green functions. 

An example illustrating the practical applications of the rules collected in this Section can 
be found in Appendix There we explicitly derive all relevant STIs for a process involving 
two gauge fields and one scalar matter field. This general analysis covers, for instance, the 
processes H — > W~^W~ , H — »• ZZ and H —>■ Z'j. Also the identities for two-point functions 
with gauge fields and scalars are discussed which will be used in our examples of Sections p.2| 
and |2]^. The drastic simplifications of that analysis within the Background Field Method 



(BFM) will be discussed in Section p73. 



2.2 Example 1: Green Functions and STIs for ^ 77 

In this Section the decomposition of the S-matrix elements in terms of IPI functions is 
described for the process if — 77 in two-loop approximation. The necessary STIs which 
relate the finite parts of the Green functions at the one- and two-loop level are discussed. 

The decomposition of the truncated, connected off-shell Green functions in terms of IPI 
functions is given at the two-loop level by the following equation: 

^ha,aM^ <12) + ^ha,aM, q2) + T%{qr + g2)G°(gi + q2)T^H\AS<li^ ^2) 
+T%\^^{q,,q,)GlM2W^ZASq2)+T^ZAM^K^^^ (2.8) 
where the tree-level propagators for the photon and the Higgs boson are given by 
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respectively. TaIap{i) is the photon self-energy at one-loop order, Thh{(1) the self-energy of 

the Higgs boson and ITij^^A^lO'i) Q^) and T^^\^^^{qi, 52) are the one- and two-loop corrections 
to the i?77 vertex, qi and q2 denote the in-going momenta of the two photons. 

In this calculation at two-loop order only QCD corrections are considered. Thus it is 
convenient to decompose the one-loop vertex corrections into a fermionic and a bosonic part 

r^H\MQi,Q2) = rgy;r(?i>?2) + rgiX(^i'^2). (2.10) 

Furthermore the two-loop terms are split into QCD and electroweak corrections: 

irgU.(^i'^2) = rgi^Jf(gi,g2)+rg^;-^(g„g2). (2.11) 

Actually, since at two-loop level only QCD corrections are considered, the terms involving 
the photon or Higgs boson self-energy in Eq. ( |2.8| ) vanish. Their contribution would be of 
the same order as the two-loop electroweak corrections to the genuine vertex. In Fig. ^ some 
sample diagrams of the remaining contributions are pictured. 




Figure 2: One- and some two-loop diagrams contribution to H —y 77. The dashed lines 
correspond to the Higgs boson, the wavy lines to the photons, the curly ones to the gluons and 
the straight ones to the quarks. 



The physical amplitude is calculated via a projection on the physical states 



-MS,, = G'g;,A.(?i,?2)6^(gi)e'^(g2) 



(2 12) 

where is the Higgs boson mass and e^(g) denotes the polarization vector of the photon 
with momentum q. 

The mass shell projection of the two-loop amplitude can be correctly performed only 
if the self-energy of the photon, r^]^^(g), satisfies the well-known transversality condition 

g'^r^^^^(g) = 0. However, in a non-symmetric regularization scheme this property is in 
general not valid any longer. It has to be reestablished as will be explained in Section ^. 

In order to obtain the complete set of one-loop counterterms, we observe that the regular- 
ized two-loop Green function r^^^^_^(gi, ^2) contains three different sub-divergences which 
require proper subtraction. Furthermore, since in general the regularization procedure breaks 
the symmetry, we are forced to introduce the following general counterterm Lagrangian 

'Ch^^^ = ^Zie QiA^qa^qi + Z2qi ^qi + Zm,Z2miqiqi + ZY,Z2YiHqiqi , (2.13) 
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where H, A^, qi and qi are the Higgs boson, the photon and the fermion fields, respectively. 
The parameters e, rrii and Fj are the gauge coupling, the masses and the Yukawa couplings 
of the fermions. In a symmetric regularization scheme the free parameters Zi, Z2, Zrm and 
Zy. are related by means of QED-WTI0. In a non-invariant regularization scheme this is not 
true any longer and these coefficients have to be fixed separately. For example, in Analytic 
Regularization which we will use for the practical computation, the renormalization 
constants are Laurent-expanded in powers of the regulators. In our case it turns out that 
it is enough to introduce only one, A [^]. Then the renormalization constants read Z^ = 
X]n>o^i"V-^" i'^ = l,2,mi,Fi). The pole parts are removed by means of the minimal 
subtraction scheme and the finite parts, Z^\ are fixed by the STIs. 

In order to fix the one-loop counterterms explicitly, we have to consider the Green func- 
tions r^^'^^^(p, p), r^^^^'^(p, and T^^l'l^'^^ (p) which arise as sub-diagrams of the two- 
loop graphs. The fermion self-energy JT^^Jj'^^'^^ (p) contains the two independent parameters 
Z2 and Zfm which are not constrained by any STL They can be fixed as usual either mini- 
mally, i.e. ^2°^ = Zj^^ = 0, or by imposing on-shell renormalization conditions 



r 



W,QCD 



0, 



d 



{1),QCD 



(2.14) 



The parameter Zy. is related to the mass renormalization constant, Z^m- In our specific 
example, where fermion mixing is absent, both parameters can be identified. Zi has to be 

{p,p) to the fermion self-energy. 



fixed in terms of the STI which relates the vertex F 



(l),QCD 



By differentiating the identity w.r.t. the photon ghost field ca and the fermion fields qi and 
g,, one immediately gets 



5ca{-P - p)5qiip)6qi{p) 



0=0 



0, 



(2.15) 



which is equivalent to the simple WTI in QED. Note that this equation is only true as 
exclusively QCD corrections are considered at two-loop order. 

After the one-loop counterterms are fixed, let us now focus on the Green function 



1,2) which is our prime interest. We again have to make sure that 



the finite parts of the counterterms are correctly fixed according to the STI. In fact, since 
the process if — 77 has no tree-level contribution there is no free overall parameter which 
fixes the finite parts by using renormalization conditions. 

According to rule 3 of the previous subsection, we consider the derivative of iS(r) = 
w.r.t. the photon ghost field ca, one photon and the Higgs boson H. As a result we 



^Z,„; is not related to Zy^ by a STI or WTI. However, in the general electroweak case two out of the 
three parameters Z^. , Zy and v, the vacuum expectation value, can be chosen independently. 
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obtain the identity which involves among others also the Green functions T ha^ {(11,(12) '■ 



<t>=o 



6cA{qi)6A^{q2)6H{p) 

+ (-icwswqi + CwT ^^^*,i{-qi)) T z^A,,H{q2,p) 

+ ^CAW;-'Hi(l2,P) {cw^Z.aM - Sw^A,A,iq2)) 

+ ^caw;-^aM ^2) {(^w'S^'z.Hip) - swTa,h{p)) 

+ ^CAG-'0i-qi)'S^GOA.Hiq2,p) + ^CAHG*-oiP, ^2)^00 A. (^2) 
+ ^caG*^OA,{P, q2)'S^G0H{p) + ^CAH*{~<ll)'^HA,H{q2,P) 

+ Tc^HH*{p, 92)^^^,(92) + Tc^h*aAp^ q2)'S^HH{p) 
= 0. (2.16) 



Actually this equation constitutes a special case of the identity (p.4| ) derived in a more 
general context. 

In the following we demonstrate how this equation simplifies for the special kind of 
corrections we are interested in. In order to disentangle consistently the QCD corrections 
from the Green functions appearing in the STI at a given order one can take the derivative of 
the latter w.r.t. the parameters of the SU{3) colour group, namely Ca and Cp. Furthermore, 
we can disentangle the contributions of the fermion loop from the contribution of the bosonic 
corrections as the coupling of ghost fields (as well as the external BRST sources W*'^, and 
H*) to the fermion lines occurs for the first time through two- loop electroweak interactions. 
Note that the Green functions ^ZpH, ^ha„, ^goa^h and ^g^h vanish at tree-, one- and 
two-loop level as they violate CP invariance. It is well-known that the CP symmetry is 
violated in the SM only through the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Thus 
CP violation manifests itself in the scalar sector starting at the three-loop order. The last 
term in ( ^.16] ) vanishes if one restricts the analysis to fermionic contributions and their QCD 



corrections. Taking these simplifications into account we finally get for the first term of the 
r.h.s. of Eq. ( [2.16 ) expanded up to two loops 



Note that the three-point function F^^^ j:^(g2,p) is absent at tree level. Further simplifica- 
tions occur through the observation that the second term on the r.h.s. of Eq. (|2.17|) does 
not give any contribution since there is no room for QCD corrections. Please note that the 
one-loop Green function r|^^^^*,3 does not contain any fermionic loop. In the same line of 

reasoning all terms except the ones in the first two lines of Eq. ( p.l6|) drop out. 

From the Lagrangian one obtains ^^^\ t.si—qi) = iswqi- This in combination with the 
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accordingly simplified remaining terms of Eq. ( |2.16| ) finally lead us to the following STI 



-<nXi^"fe,p) = 0, . = 1,2. (2.18) 

This identity must be fulfilled at one- and two-loop order for the specific corrections we are 
interested in. 

In Section ^ the breaking terms for these identities will be provided. Furthermore we 
will compute the amplitude in the analytical regularization scheme and we will show how 
the algebraic renormalization works for this two-loop example. 



2.3 Background gauge fixing 



As is well known the BFM |Tl]] allows to derive the S-matrix elements in terms of Green 
functions with external background fields with the exception of fermion fields. The main 
simplification in the BFM results from the fact that the theory with background fields 
possesses two different invariances: the ERST symmetry, which involves quantum fields and 
ghosts, and the background gauge invariance. The latter provides several simplifications 
in the computation of physical amplitudes and in the renormalization procedure due to its 
linearity. The first systematic application of BFM in the SM for invariant regularizations at 



the one- loop level was presented in p5|. Our considerations regarding BFM, however, also 



apply in the case of noninvariant regularizations and also beyond the one-loop level. We 
focus on the practical aspects of the BFM which are relevant for higher-loop calculations. 
Further details on the theoretical advantages of the BFM can be found in p5|, ^ pl|] . 

The main difference between the STIs ( p.6| ) and the WTIs for the background gauge 
invariance is due to the linearity of the latter. Linearity means that the WTIs are linear in 
the functional F and therefore they relate Green functions of the same orders while for the 
STIs there is an interplay between higher and lower order radiative corrections. 

To renormalize properly the SM in the background gauge, one needs to implement the 
equations of motion for the background fields at the quantum level. The most efficient way 
to this end is to extend the BRST symmetry to the background fields 

= 0, 

sn^ = 0, (2.19) 
sn^ = 0, 

where and fi^ are (classical) vector fields with the same quantum numbers as the 

gauge bosons W, Z and GJJ, but ghost charge —1 (like an anti-ghost field), fi^,^'' and 
are scalar fields with ghost number —1; in the following we will denote by Q the complete 
set of these fields. 

In the following we will denote with F' the effective action which depend on the Q (and 
correspondent ly the its Slavnov- Taylor operator) and with T = T'\^^q, i.e. by setting Q 
to zero. 







s^l = 


0, 


5(7° = 








s^t = 


= 0, 


sG^ = 






= ^;, 




0, 


sH = 
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Therefore one has to modify correspondingly the STI 



s'iT'M = s{rM + 1 d'xlnl 



'6T' 6T'' 



Cw 



.6Z„, 5Z, 



- Sw 



' 6T' 6T' 




.6A„, 



AGO SG\ 



5H SH 



(2.20) 



Here iS(r')[0] are the STIs given in Eq. ( |2.6| ). Notice that also the scalar fields G"^, G° and 
H are paired with their own background fields, G^, G^ and H in order to extend the 't Hooft 
gauge fixing to a background gauge invariant one described in Appendix 0. 

To study how the effective action F' depends on W^^, for instance, one has to derive the 
STI ( p.2q ) with respect to the fields After setting fl^ = one obtains 



6T 



6T 

6Wf 



T 



(2.21) 



where St is the linearized (conventional) Slavnov- Taylor operator of Eq. (A. 8). These equa- 
tions describe the relations between the quantum and the background fields and they sup- 
plement the STI and the WTI: 



5(r) = o, W(A)(r) = o 



(2.22) 



where W(a) is the WTI operator of the background gauge invariance (cf. Eq. ( |2.24|) ). 

The space of counterterms and of possible breaking terms to the STI is enlarged by 
those monomials which contain the background fields and the fields Q in addition to the 
conventional fields and anti-fields 0*. This requires a new analysis. 

In the calculation of the necessary counterterms, there are two main approaches. They 
only differ by the ordering in which the three equations in (2.21) and (2.22) are used. 

1. We first use the WTI for the background gauge invariance and a subset of the conven- 



tional STIs (|2.6| ). The only missing parts are counterterms which relate the two-point 
functions of the background field to the two-point functions of the quantum fields and 
for the two-point functions r^*Q, where 0* is a generic anti-field. To fix these last 
counterterms one has to use the extended STI ( |2.20| ). 

Thus, to fix for example the counterterms for the two-point function for the quantum 
field W^, one has to derive the STI (|2.20|) with respect to WJ" and with respect 

to n^,w^: 



T 



■ip) 



w+w- 



T 
T 



w+w- 



(p) + Tn+w'.-ip)Tw+w-{p) 



w+w 



ip) + Tn+w*,-ip)^ 



w+w- 



(2.23) 
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These equations fix completely the two-point functions r^+^_ (p) , r^y+^_ (p) in terms 
of the quantum one r^+^~(p). Finally, in order to fix the counterterms for the two- 
point functions r,^*^ (a complete discussion has been given in |^), one has to consider 
the derivative of (|2.21| ) with respect to the anti-field 0* and one ghost c. This conclude 



the algebraic renormalization program in the case of the BFM. 

2. The second alternative approach exploits completely the use of the extended STI, 
namely equations of the type (|2.21|) : Thus, one first computes all possible counterterms 
to restore the WTI, then one fixes the remaining counterterms by considering the 



functional derivative of the extended STI (|2.20| ) with respect to fi, the background 
fields (f) and the quantum fields 0. However, as in the former approach this cannot 
exhaust completely the algebraic renormalization program: One shows that one still 
needs a reduced set of STIs in addition which can be derived from the conventional 
STI (|2.6| ). In particular, besides the Eqs. (|2.21|) , one needs the STI to fix the anti-field 



part of the action. This guarantees that the BRST transformation are preserved to all 
orders. 

In a practical analysis of a physical process where one has not to compute the same Green 
functions with the quantum fields replaced by the corresponding background fields nor vice 
versa, the first approach is favourable — as in the two phenomenological examples discussed 
in this paper. One can disregard all equations of the type ( p.21|) . This simplifies the analysis 
significantly. A phenomenological example using the second approach will be discussed in 
a forthcoming publication. In Section |3.2| we discuss the analysis of the breaking terms in 
both methods within the BFM. 

In order to single out the relevant WTI, we have to take into account the rules stated 
in Section However, in the case of the background fields, the role of the ghost particles 
is played by the parameter of the infinitesimal background gauge transformations (see Ap- 
pendix To each generator of the gauge group 5'?7c(3) x SUj{2) x ?7y (1) we consider the 
corresponding local infinitesimal parameters. They are denoted by \a{x), \z{x) and A±(x) 
for the electroweak part and Xa{x) for the QCD sector. Thus, the functional WTI for the 
effective action F' reads: 

/rip/ „ 
jir-^= d'xK{x)w'^{x){r) = , (2.24) 

<p 0(p{X) a=A,Z±,a-' 

where the variations 5x(f){x) are explicitly given in Appendix ^ (see Eqs. (|A.4| )-( [A77| )). The 
sum runs over all possible fields and anti-fields. W^(A)(>V4(a;)) is called Ward-Takahashi 
operator and acts on the functional ir'[0]. An explicit expression is given in 



Concerning the rules of Section only two slight modifications of the rules ^ and ^ are 
necessary: 

3'. One has to differentiate the general WTI ( |2.24| ) w.r.t. the infinitesimal parameters 
Xy in order to get constraints on the Green functions involving the corresponding 
background gauge fields V . 
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4'. To derive constraints on Green functions involving one ghost and one anti- field plus 
other quantum fields one either can derive a corresponding STI with rule ^ (i.e. differen- 
tiate Eq. ( p.6| ) w.r.t. two ghost fields) or one can derive a linear WTI by differentiating 
w.r.t. one ghost field and one infinitesimal parameter A^. We prefer to use the second 
version since linear WTIs are simpler to handle within our specific subtraction method 



(cf. Section |3]^). Note that this choice implies some assumptions on the wave func- 
tion renormalization for multiplets of fields as will be explained in more detail in the 
example on 6 — > 57 (cf. Section 

In connection to these modifications a remark is in order. If we consider a Green function 
with one gauge field, one ghost field and one anti- field (e.g., ^caW*-+w-) have to differ- 
entiate the WTI ( p.24| ) w.r.t. the ghost field, the anti-field and the infinitesimal parameter 
(which in the example is A_) associated to the gauge field. This provides an identity which 
fixes the considered Green function. In the case that no gauge field is involved (e.g., Tcj^q*q') 
one has to consider the background variation of the ghost field, of the anti-fields and of the 
quantum field (which in the example is q'). Thus the WTI has to be differentiated w.r.t. 
^aCa, q* and q', and ca, S\q* and q', and ca, q* and 6xq'. Some of the resulting WTIs coin- 
cideQ. One has to select the independent ones, but this can be easily done by inspection of 
the WTIs themselves. 

At this point let us consider the example discussed in Appendix |B] the amplitude involving 
two gauge fields V/^ and one scalar field $ in the context of the BFM. Thus we are able to 
compare the two approaches and to underline the differences. 

In the framework of the BFM the two gauge fields, Vi and V2, and the scalar field $ 
are replaced by their counterparts V"/*, V2 and respectively. As in the conventional gauge 
fixing the amplitude is built up by irreducible Green functions which in this case read JTyt^yi, 

= 1,2), and Ty^y^^. In the following we will denote irreducible Green functions 

where only external background fields are involved as background Green functions[]. 

Let us in a first step consider the two-point function r^i^i(p). We get the following 

identities using ( |2.24|) in combination with Eqs. (|A.4| )- ([A77|) : 



5'W(A)(r) 



5Xy{-p)5VS{p) 



ip^Ty,y3{p) + Y.M,^^,T^,y3{p) = 0. (2.25) 



The sum runs over all Goldstone fields and with masses M-j. = i^M^y, M^ qo = 
—Mz and zero for all the other combinations. In the following the summation sign will 
be omitted. A comparison with the corresponding identity in the conventional formalism, 
Eq. (|B.1|) , shows that Eq. (|2.25| ) is linear in the Green functions. However, it requires the 



renormalization of the mixed two-point functions W^iyoip) which can be studied with the 
® This is a consequence of the consistency conditions to be discussed in Section 

^Notice that in the foUowing equations where we consider single components of WTIs or STIs, or even 
some specific Green functions obtained from IT' we can avoid the prime. 
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help of rule 3': 



5'>V(A)(r) 



(2.26) 



<j>=0 



Let us now come to the three-point function Tyi^y^^{pi,p2). From Eq. ( p.24| ) we get 

5'W(A)(r) 



SXV{-Pi-P2)SVj{p,)6<l>{p2) 



</>=0 



i (Pl +P2)Ty^,^.^(pi,P2) + Mi^^,T^,y,^{pi,P2) + fijk^vk^{p2) + ti,<I.$'r^,^,(pi) 

= 0, (2.27) 

where fijk and represent the structure constants, respectively, the generators of the 

gauge group in the representation for scalar fields. We refrain from listing them explicitly. 
In the above identity the two-point functions are already known and only the function 
ri,^ji (pi,P2) is new. It is fixed by the WTI 



5'W(A)(r) 



6X' 



vK-Pi -P2)S<l>{pi)6<l>'{p2 
i{p 





1 +P2)Ty.|,,|,(pi,P2 



0=0 



i>'(pi,P2) + ii,<i.'<i."ir|,„|,(pi) + 1 



(2.28) 



where again the sum over $" takes the values and G^. 

The four equations (|2.25|) , ( p.26|) , ( |2.27|) and (|2.28| ) already form the complete set of 
identities needed for the computation of the amplitude. In fact, all identities are linear in F 
and therefore they keep the same form to all orders. This also implies that the coefficients 
fijk,ti^<s>q,' and M( <j, are not renormalized. Their renormalization is fixed from the renormal- 
ization conditions. Furthermore no Green function involving ghosts or anti-fields occur. Let 
us mention that instead of the four identities derived above roughly ten mostly non-linear 
STIs have to be analyzed in the conventional gauge fixing. Thus, in this case the BFM is 
obviously superior as compared to the conventional gauge fixing. 

At this point a practical remark is in order: the obvious advantage of the BFM due 
the linearity of the WTIs is only valid at the highest order of the computation. In lower 
orders, i.e. in sub-diagrams, also quantum field Green functions are involved which makes 
it necessary to use all three types of identities, namely (|2.20| ), (|2.24| ) and ( p.21| ), in general. 
Thus, in specific examples, the BFM may introduce more complications in the sub-diagrams 
in comparison with a conventional gauge fixing such that the advantages of the BFM at the 
highest loop level could get partly compensated. However, in the example of the two-loop 
corrections to 6 ^ S7, to be discussed in the next section, we will show that the analysis 
of the sub-diagrams within the BFM is still favourable compared with the analysis within a 
conventional gauge. 

Let us finish this section with two practical remarks about the gauge fixing and renor- 
malization conditions within the BFM. 
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In order to evaluate the S-matrix elements in the BFM a gauge fixing has to be chosen for 
the background gauge fields. However, this choice is completely independent from the gauge 
fixing used for the internal gauge fields. This allows for a more convenient choice oriented 
on the physical process. For instance, the BFM Green functions with external unphysical 
scalar bosons and G^), with external ghost fields as well as longitudinal gauge bosons 
can be neglected in the decomposition of S-matrix elements in terms of IPI parts. This can 



be achieved via the use of the unitary gauge fixing for the BFM propagators (see pSj for 
explicit examples). The procedure can be implemented both for the electroweak and the 
QCD sector of the SM. 

Besides the symmetries of the BFM one has to impose some renormalization conditions 
in order to unambiguously fix the finite parts of Green functions. It appears very useful 



to implement them in terms of background Green functions . The relation between the 
renormalization conditions for the background Green functions and those for the quantum 
fields are considered in [ETI . 



2.4 Example 2: Green Functions, STIs and WTIs for the 6 ^ 57 

In this section, we briefly describe the decomposition of the S-matrix elements for the process 
6 ^ 57 in terms of IPI functions at two-loop approximation. Simplifications concerning the 
renormalization procedure are discussed in the context of the BFM. We explicitly derive all 
WTIs and STIs constraining the counterterms at the one- and two-loop level for this specific 
process . 

The decomposition of the (truncated) connected BFM Green functions in terms of IPI 
functions can be split into two contributions. The first one, G^f'l^^^^^ipb.Ps)^ contains 

A^bs 

the flavour changing neutral current (FCNC) vertex corrections whereas the second one, 
G^f^'^f'{ph,Ps), the FCNC self-energies: 



Gfjp„p,) = Gf:f^^{p,,p,)+Gf-{p,,p,). (2.29) 



Remember that for the background fields we have chosen to use the unitary gauge in order 
to avoid external unphysical particles. Then the two contributions are given by: 



Au,bs \fo^fS) ^, AnA„ AaAv A^hs bb bb A^bs Anbs '^^ 



+ r^G^^rZ, + r^^Gg'r"' . (2.31) 

We recall that G\f denotes the tree-level propagators and T the irreducible Green functions. 
After projection on the physical states, the contributions from the 7 — Z and 7 — if mixings 
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vanish because of the WTI 



<5'W(A)(r) 



6Xa{-p)SH{p) 



(2.32) 



0=0 



where an analogous equation holds for F^^^^^. Thus the second line in Eq. (|2.3CI|) and the 
second and third lines of Eq. ( p.31|) drop out from the amplitudes. 

Let us in a first step consider the two-loop contribution to the 6 — 57 amplitude and 
derive the corresponding WTI. According to our rules we have to replace the photon field, 
A^, by the corresponding infinitesimal parameter ol0 Uq{1), Xa- Furthermore we have to 
take the derivative of Eq. ( |2.24| ) w.r.t. Xa, b and s. Using the formulae ( [A.4| )-( [A77D we 
obtain 



<5'>V(A)(r) 



(2) 



4>=o 



SXa{-P.s - Pb)Ss{ps)6b{pb] 

KPs+PbrTfjps,pb) + teQ, (rfJipb) - rg)(-p, 



n 




(2.33) 



where Qa is the charge of the down-type quarks and pb and Ps are the in-going momenta of 
the quark lines. This identity (for a one- loop analysis see also |2^) can be used to fix 
the overall counterterms defined through 



££f/ = zg;iAS7^PL6 + ^PLb + z'i^Axsj^Pnb + Z^'h ^Pnb + h.c. , (2.34) 

where the Z factors, in general, contain finite and divergent contributions. Clearly the same 
Lagrangian also holds at one-loop order. Pl/r = (l=F75)/2 are the projectors on the left- and 
right-handed components. Note that for an invariant regularization scheme no counterterm 
at all is needed for the Green function ]r'''^\ (pb,Ps)- However, if the regularization scheme 
breaks the identity ( p.33| ), it can be restored with the help of non-invariant counterterms 
in (PI). 

Let us now have a closer look to the sub-divergences. From the topological structure of 
the two-loop diagrams with only external background gauge fields, which are contained in 
r*^^'* (cf. Fig. |3|), it is evident that the IPI three- and four-point Green functions with 
external quantum, respectively, ghost fields do not appear at one-loop order: 

-[p(l) Tp(l) Tp{l) (2 QK] 

Here are the vector quantum fields and c" are ghost fields. Also the Green functions where 
the vector fields are replaced by scalar fields do not contribute. Actually the renormalization 
of sub- divergences with more then two quantum fields enter the calculation only at three- 
loop order. Thus, we only have to consider Green functions which are either sub-diagrams of 

®Here the index Q reminds that the abehan group of QED is meant. 
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Figure 3: One-loop diagrams contributing to 6 ^ 57. In order to obtain the two-loop diagrams 
one-loop corrections to each vertex and internal propagator have to be considered. The two 
external fermion lines represent the b, respectively, the s quark and the external wiggled line the 
background photon, A^. 



bs diagrams occurring in (|2.30|) and ( p.31| ) respectively. Concerning the three-point 



diagrams the following functions have to be taken into account: 

Tf^-JPs,Pb), Tf^.^Jp2,pi), Tf^^^^^.{p+,p_), Tf^^^^_{p+,p.), 

(2.36) 

^aIg+w-(p+^P-^^ ^aIw+g-(p+^P-^^ rgi-^^^(p2,Pi). 
The two-point functions 

{Ps,Pb) , ^'iJq, {P2, Pi) , ^w+w- (P+^P-) ' 



^G+G-iP+^P-) ^ ^G+W-^P+^P-^ ^ ^W+G-^P+^P- 



(2.37) 



appear as self-energies in the two-loop graphs. In Eqs. ( ^.3(j| ) and ( ^.37] ) qi and q2 are two 
generic quark fields. Notice that the Green functions F^^^^+j^- and T^^^^+q. are absent at 
tree level. This is a consequence of the choice for the gauge fixing and the background gauge 
invariance. At one-loop level, however, contributions may appear as soon as a non-invariant 
regularization scheme is used. 

By inspection of Eqs. ( p.36| ) and ( |2.37| ), we see that the same Green functions with the 
quantum fields replaced by the corresponding background fields never occur. This implies 
that the STI identities ( p^.2lD relating background and quantum fields are irrelevant for the 
practical analysis of this specific process and we can restrict ourselves to the conventional 
STI and the WTI. In Section |3.2| we present a general discussion on this point. 

Thus, let us discuss the STIs and WTIs which constrain the relevant sub-diagrams. Using 
our rules for the BFM we are able to derive the complete set of identities. For the three-point 
functions involving background photon field WTIs are used whereas only a reduced set of 
STIs for three-point functions and two-point functions are indeed necessary. 

We start with the WTIs. From ( p. 241 ) one gets: 



53n;(A)(r)(i) 



SXa{-Pi - P2)Sq2iP2)Sqi{pi 



0=0 
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(2.38) 



which is the analogue equation to ( [^.331 ). Here, however, qi and q2 refer to any type of 
quark fields. As already mentioned above, an advantage of the BFM is the linearity of the 
identities w.r.t. F. This allows to disentangle easily the fermionic corrections from the 
bosonic ones. In fact, in the same way as for the H 77, we can select the independent 
contributions and introduce the corresponding counterterms. Notice furthermore that no 
ghosts are involved. In the case of conventional gauge fixing the identity ( p. 38 ) would be 
replaced by a STI obtained by differentiating w.r.t. ca, qi and g2- Already at one-loop order 
this identity would require the computation of Green functions involving ghost particles and 
off-shell quark fields. 

The other Green functions of Eq. ( p.36|) involving background fields are constrained by 
the following WTI 



(1) 



5\a{-P+~P-)6W+{p+)6W-{p. 



<j>=0 



I {p+ + p^Y ^%w+w^ ip+,p-) - «e ( rJ^^V,,,- {p- 



0, 



(53>V(A)(r)(i) 



8\a{-P+-P-)6G+{p+)6G-{p 

AaG+G 



(t>=0 



0, 



5^>V(A)(r)(i 



8\A{-p^-p^)m^{p^)bG-{p^ 



0=0 



0, 



-P^ 



(2.39) 



(2.40) 



(2.41) 



and their hermitian counterparts. 

According to our third rule of Section for the Green functions r^^^+^^-(p+,p_), 

^i^iy+G-lP+'P-). A^G+w-iv+^P-) and T ^^^^^.{p^.pS) of the above equations one has 
to derive new STIs. For instance, one gets 



(525(r) 



M^(-p+-p_)5c+(p+)5W^-(p_) 

r 

= 0, 



?i=0 

A^c+wt^-^V+.P~)^w^w-^V-) + ri^c+G*.-(P+'P-)rG+iy-(P-) 

p+)^A^w'+H/-(p+>p-) + rc+G*>-(-p+)ir^^G+i^-(p+'P- 



(2.42) 
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where the new Green functions ^+^*~{p^,p^) and r^^^+g,.,_(p+,p_) emerge. Clearly, 
also the STIs ( |2.42| ) can be spoiled by the radiative corrections and therefore it must be 
restored by suitable counterterms. However, it turns out that due to the consistency con- 
ditions (see Section such STIs deliver no independent constraints. They are auto- 
matically preserved if the STIs for the two-point functions (discussed below ( p.48| )), the 
WTIs (p.39|) - (p.41|) , and the WTIs involving the new Green functions r^^^+^.,-(p+,p_) 
and ^+^,,_(p+,p_) have been restored. These new WTIs are obtained by differentiating 
the WTI ( 12:2^ ) w.r.t. Aa,c+(j9+) and W*^~{p^) and w.r.t. Aa,c+(p+) and G*'-(p_): 



5'W(A)(r) 



(1) 



(/>=0 



0, (2.43) 



53>V(A)(r)(i) 



5\a{-P+ - P-)5c+{p+)5G*^-{p^) 



0=0 

(1) 



I (p+ +p_)'^r«^^^,,_(p+,p_) - le {t%,AP-) - ^'•clG*^A-P^ 
= 0. (2.44) 

These WTIs do not involve any new Green functions since the two-point functions r|^^''^*,_ 
(1) 



and r|:+g.^_ are already fixed by means of the Faddeev-Popov equations (|B.2|) as discussed 



in Appendix |^. Notice that only the first equation of ( p.43| ) can be broken by the radiative 
corrections since the second one involves (by Lorentz invariance) only finite quantities. 

To restore the identities of the sub-diagrams, we need a complete set of counterterms. It 
is convenient to divide them into three different sets 

^(1) _ Al),WTI Al),STI Al),INV (r, 

organized in such a way that 

w,«(/dMir) = w,„(/dMir) = •). 

5o(/d\4ir) = 0. (2.46) 

This triangular organization of the counterterms ensures that it is possible to restore the 
WTIs by only fixing the coefficients of C^^}^^'^^ . The counterterms C^f^^^^^ and 

are invariant under the background gauge transformations. is necessary to restore 

the STIs and with the help of the renormalization conditions can be fulfilled. In 

Section |3]^ and in Appendix we will prove with the help of the consistency conditions that 
this procedure is always possible provided no anomalies occur. 
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The complete list of counterterms for three-point functions needed to restore the WTIs 
is given by 

+ {e.^xaA^W^'^^V'^W-''' - e,,xaA^W~'''V''W+''') 

+ Z^'^AxG+V^G- + z['^AxG-V^G+ + Z^^AxG+Wx 
+ zil^AxG-W+ + [zil^AxsPLi'b + zi^AxsPR-f^b+ 

+ zil^Axq2PLl\i + zSAxq2PRl\i + h.c.) 

+ Z^^Axc-^W*x'~ + zSAxc-W*x'^ , (2.47) 

where V is the covariant derivative w.r.t. Uq{1). These counterterms correspond to the vari- 
ous Green functions occurring in the WTIs (P^)-(|CT|) and in the WTI {^^. In Eq. (^^) 
they are partially written in covariant form w.r.t. Uq{1) and also include counterterms to 
the WTIs which, however, are not relevant for our specific process under consideration. They 
would be needed, e.g., for the calculation of four-point functions like r^^^_^^+^y-. We note 
that this is the preferable basis for counterterms in order to analyze the renormalization of 
the whole model |21[ . 

Besides the WTIs, we also have to take into account the following STIs. The general 
form for the two-point functions is also discussed in Appendix ^ In Eqs. ( [B.l| ) and ( [B.2| ) 
they are given for a generic gauge field, V^, and a generic scalar field Note, that the 
Faddeev-Popov equation has to be used in order to obtain this simple form 



6c+{-p)6W-{p) 



6c+{-p)5G-{p) 



<j>=0 



r^+^..-(p)r^+e.(p) + r,+G..-(p)rG+G-b) = o, (2.48) 

which at one-loop order become 

r(^4,,4p)(2M^^7.M) + *p^r(;V^^_(p) + rf4^^ = o, 

T'^^^^..{p){-tpMw) + ^P.r^^+a-(P) + rl^G'^^^^^ = 0- 



(2.49) 

Notice that the ghost fields do not couple directly to fermions. Hence, at one-loop level these 
identities can be separated into two sets. In fact, by decomposing F^^^ = ]rW'f^^"^ -\-]rW''^^^* ^ 
where ]rWJerm contains only diagrams with virtual fermions and r*^^^''"^'** the remaining 
ones, the fermionic contributions satisfy the following simplified identities 



'^^fr(p)+'M^^c4 
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^P^r(iV;:-(p) + zMv.rgi'^T-(p) = 0. (2.50) 



These identities have to be considered if the computation is done in the 't Hooft-Veltman 
scheme. 

In general, the two STIs in (|2.48|) can be restored by counterterms which are invariant 
under the background gauge transformations (That this is possible is shown in Appendix |C[). 
The corresponding Lagrange density reads: 

+ Z^iv^'G-W^ + Z^^V^G^V^G- + zg^G+G" 

+ z^^lhPL W^q + z^k'PL W^s + h.c. . (2.51) 

Note that the last two terms are background gauge invariant because the W boson transforms 
as an isovector of SU{2). 

Also the following STI for the three-point functions have to be considered 



535(r) 



4>=o 



^c+W:--iP'l + Pb)'^W+gbiP<J^Pb) + ^c+G*:-{Pq + PbWG+qb{Pq,Pb) 



Tgg,{-Pg)T^+gUt,{Pg,Pb) " T ^+ {Pg , Pb)T g, b{p 



'b 



= 0, (2.52) 

where the sum over the quark fields q' is understood. Here q and q' are generic quark fields 
and g* and g'* are the corresponding BRST sources. An analogous equation where c"'",g 
and b is replaced by c~ ,s and q has to be taken into account. We only need the one-loop 
expansion of these identities which reads for ( |2.52|) 

= {Pg + Pb) iiYgbluPh) +iij>q+ Pb)^ rj^Vg-fe(P9, Pb) 

+ r^i.)^.,. {pg + Pb) {i^ [PLnig - Pnmb]^ + t MH.rgl-,(p„ pb) 

- z(- ^g - mg)T%,,{pg,pb) " T ^ ( -J), ) (-^X^.^^Pl) 

- {tVgg^Pi,) r'-^lipb) - T%,.{pg,pbWb ~ nib) • (2.53) 

Here Vgg' are the CKM matrix elements. For convenience the prefactor e/{swV^) has been 
omitted. 

As stated above, the study of the one-loop approximation disentangles the different con- 
tributions coming from fermionic and bosonic radiative corrections. Unfortunately in the 
case of Eq. ( |2.53| ) it is very hard to disentangle the fermionic contributions because of the 
presence of external fermionic fields. In addition we also have to note that the Green func- 
tions with external ghost and anti-fields — in contrast to the case of the STIs for the W boson 
two-point functions — contain fermion loops already at one-loop order. This is because there 
are vertices involving the anti-fields q'* and q'*, ghosts and fermions (see |18[ and 

Appendix ^ for the Feynman rules). 
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If we consider on-shell b quarks some Green functions vanish and the identity ( p.52| ) is sim- 
phfied. This can be heavily exploited in the algebraic one-loop analysis of 6 ^ 57 How- 
ever, we have to remember that since these one-loop corrections appear as sub- divergences 
for two-loop amplitudes these simplifications do not apply here. 

In Eq. ( |2.53| ) Green functions with fermionic BRST sources are involved, like, e.g., 
^^^+-itf^{Pq,Pb) and their hermitian counterparts. Thus, according to rule ^ (cf. Section pTT| ) 
one has to consider STIs for them. However, in the case of the BFM, according to rule 4', 
we are left with the following linear WTIs: 



<5"W(A)(r) 



6XA{0)Sc+{-pq - pb)Sq*{pg)6b{pb) 



0=0 



- ieTc+q*b{Pq,Pb) + ieQqTq*c+b{-Pq " Pb,Pb) " IcQbTbc+q* {-Pq - Pb, Pq) 



^'W(A)(r) 



5XA{Gi)5c+{-pq - pb)5qi{pq)5h*{pb 



<j}=0 



- ieTc+g^b*{pq,Pb) + ieQqTg^c+b*{-Pq - Pb,Pb) - ieQbFb*c+qA-Pq - Pb,Pq) 

= 0. 

The corresponding identity where, e.g., Aa(0) is replaced by Xz{0) reads 



(2.54) 



5'>V(A)(r) 



6Xz{0)Sc+{-pq-pb)6qt{pq)6b{pb) 



<t>=o 



- Mz^ff^c+q^bi-Pg - Pb,Pq,Pb) 



+ ie—Tc+q*b{Pq, Pb) + ieQ^ Tq*c+b{-pg - pb, pb) 
Sw " 

0. 



ieQf Tbc+q* {-Pq - Pb, Pq) 



(2.55) 



Here Qf/^ = {Qu/dSw I cw T ^1 {'^swCw))- There are two other pairs of equations where 



Aa(0) and 
respectively. 



-Pq — Pb) are replaced by A+(0) and CA{—Pg — Pb) or A+(0) and cz{—pq —Pb), 



Finally we have to fix the invariant counterterms |]I6| , |T^, [T^ C 



,{1),INV 



in order to fulfill 

the renormalization conditions. It contains Zw, the wave function renormalization for the W 
boson, its mass M^/, the wave function renormalization for the Goldstone boson, Zq, and the 
corresponding mass Mq (which coincides with the product of the gauge parameter and 
M^). Furthermore we have to fix the renormalization conditions for the fermions, namely 
the masses rriq- inside of loops^ and and which are needed for the computations of 
the on-shell amplitudes. Also the CKM elements and the couplings aqEo and Gp have to 
be fixed. 

For the CKM matrix there are two possible choices which can be adopted [^: (i) the 



use of the MS scheme where only the poles are subtracted [3C, 31 1 or (ii) the definition 



^As already mentioned in Section |2.3| rule 4', the fermion wave function renormalization is fixed by 
WTI ( 2.54 ). For further discussion see ||2l|. 
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given in which rehes on subtractions at zero momentum. For the general analysis of 
renormalization conditions in the background field gauge we refer the reader to I^Tj. Our 
specific choices in the case of b —>■ s'j will be discussed in Section |^. There we will see that 
some of the identities are automatically preserved by a conscious choice of regularization. 
This will provide great simplifications. 



3 Renormalization of the identities 



In Section P . 1| we give a brief review to the Quantum Action Principle (QAP) from a practi- 



cal point of view and introduce some necessary notation. We present the principle algebraic 



procedure necessary to remove the breaking terms. In Section |3.2| we briefiy discuss the im- 
portant practical consequences of the consistency conditions. Finally in the last subsection, 
we propose our strategy which provides the possibility to remove the breaking terms in an 
efficient way. 

3.1 The Quantum Action Principle and the algebraic method 

The QAP is the fundamental theorem of renormalization theory. It guarantees the locality 
of the counterterms, and as a consequence, the polynomial character of the renormalization 
procedure. The QAP also implies that all breaking terms of the STIs and WTIs are local 
and that they can be fully characterized in terms of classical fields^. 

Formally, the QAP states that within a specific renormalization framework derivatives of 
a IPI generating functional F w.r.t. a parameter^ of the theory P], /z, or w.r.t. a field 
are local insertions A in the IPI Green functions 



dfi ^ 6(f){x 



^M-r, ^^=A^(x)-r. (3.1) 



The explicit meaning of the r.h.s. is the following: In analogy to ( |2.5| ) we consider an 
(extended) action at lowest order 

To = 1 d^a; (^Cinv + E '^I ^'^^ + ^gf + Uu + T. Pr^''^ ■ (3-2) 

where the sum runs over all the possible local insertion A^. Then one has 
A"^ ■ r = (3.3) 

E j— / n d^Pi J2Pk MPi) ■ ■ ■ (l^n-l{Pn-l)pr{Pn)^''^l.^^_,p^{Pu ■ ■ ■ , Pn) , 

n=0 \j=0 / \ k / 



^"^Within this subsection, we will skip the prime on the effective action and on the functional identities 
since the following results are clearly not restricted to the BFM 

Here we mean all the parameters of the renormalized theory: masses, couplings, vacuum expectation 
values, gauge fixing parameters, renormalization scales, infra-red (IR) regulators, . . . (see [Esl pj). 



25 



with 

501 (pi) . . . 50„-l(p„-l)Pr(Pn) 



= n^i.0„_.p.(Pl,---,Pn). (3.4) 

</>=0 

Thus, A*" ■ r generates the IPI Green functions with an insertion of an integrated or local 
composite operator 

(T(0i...0„A^'-))''''- (3.5) 

It can be decomposed into a basis of integrated monomials of fields and their derivative 
with the same quantum numbers as the l.h.s. of Therefore the r.h.s. of (|3.1|) can be 

decomposed into the classical insertion and their radiative corrections: 

A^-T = A^ + OihA^), A^ix)-T = A^ix) + OihA^ix)). (3.6) 

In the case of STIs, the QAP implies that the (subtracted) Green functions F, computed 
within a given scheme, fulfill them up to local insertions Asti in the IPI Green functions: 

5(r) = Asti + 0{hAsTi) ■ (3.7) 

Here Asti is an integrated, Lorentz invariant polynomial (of the fields and their derivatives) 
with ultra-violet (UV) degree < 4 and IR degree > 3 (assuming four space-time dimensions). 

Although Eqs. (|3.1|) and (p.7|) apply to any renormalization scheme, the coefficients of the 
various As depend on the particular scheme adopted. In fact, the definitions of A^, A^{x) 
and Asti rely on specific conventions for composite operators. Thus a renormalization 
description for the composite operators is necessary. Here one uses the concept of Normal 



Product Operators (NPO) introduced by Zimmermann or the conventional counterterm 
technique which is preferable from the practical point of view. 

Once the breaking terms Asti are given we can discuss the main objective of the algebraic 
method |T^. This essentially entails in a prescription to restore the identities by suitable 
local non-invariant counterterms, P*^^, such that at ra*^ order one has 



n-l 



5(r)(") = 5o(r(")) + 5](r(-'), r^"-^)) - 5o(r^^'(")) = oihAsTi) , (3.8) 

where the decomposition given in Eq. ( [A.8| ) has been used. Notice that the Green functions 
r(^) with j < n are already fixed and only P*^") has to be adjusted by the local counterterms 
^^.(n)^ 'jj^us^ in practice the problem amounts to solve the algebraic equations 

5o(P^^'(")) = A5r7, (3.9) 
where Sq is given in Eq. (^]S). This equation turns out to be solvable in absence of anoma- 



lies [|], |T0| where only the consistency conditions have to be used (cf. Section |3^). Moreover, 
due to a non-trivial kernel of the operator So (i.e. the space of invariant counterterms), one 
is allowed to impose renormalization conditions tuning the free parameters of the model. 

This principal algebraic procedure does not automatically lead to a practical advice 
for higher-loop calculations. As already mentioned in the Introduction, regardless which 



26 



regularization scheme one uses, the calculation of Asti in ( P-9| ) is quite tedious and gets even 
more complicated at higher orders. In general, one has to calculate all Green functions which 
occur in the complete set of STIs. Inserting them in the STIs one fixes Asti- The additional 
computations necessary in the conventional algebraic method can be slightly reduced: instead 
of calculating all Green functions which occur in the full set of the STIs, one can simplify 
the problem by computing the Green functions in special points, namely for zero momentum 
p = 0, for on-shell momentum or for large external momenta. As a consequence the breaking 
terms, Asti, are simply related to Green functions evaluated in these special points. Clearly, 
if on-shell renormalization conditions are used in the calculation, the on-shell method is 
definitely superior to the zero-momentum subtraction. In the infinite-momentum scheme 



one can take advantage of Weinberg's theorem (see also Section |37^ ). 

Despite of these simplifications, still all Green functions involved in the STI have to be 
taken into account for the computation of the breaking terms in (|3.9|) . In Section |3^ we will 
present our strategy which drastically reduces the additional work as will also be shown in 
the examples of Section H and ||. 

3.2 Consistency and renormalization conditions 

One of the main tools of algebraic renormalization is provided by the algebraic relations be- 
tween the functional operators S'r' of the STIs (see Eq.( p.20|) and the definition ( |A.8|) ) , 
of the WTIs (given in Eq. ( |2.24| )) and of the other supplementary identities like the Faddeev- 
Popov equations (see 0, 0). Beyond their relevance in the theoretical framework the 
consistency conditions turn out to be important for practical applications. 
The operators S'jp' and W(\-) form an algebra 

S'l, = if S'{T') = 0, 
5V(%)(r))-W(\)(5'(r)) = 0, 

>V(\)(W(%(r))-W('^)(W(\)(r)) = W(',^^)(r), (s.io) 

which, applied to the breaking terms A'^ and A'^y(A) of the STI, respectively, WTI 

S'iT') = A's + 0{hA's), 
W(\)(r') = A'^(A) + O(^AV) , (3.11) 

leads to the so-called consistency conditions 

S',{A's) = 0, (3.12) 

5i(A'^(A))-W(\)(A^) = 0, (3.13) 

>V(\)(A'^(/3))-W(%(AV(A)) = AV(AA/3), (3.14) 

where (A A /5)" = /"'"^A''/?'^. These kind of equations are called Wess-Zumino consistency 
condition. 

The consistency conditions have very important practical consequences: 
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1. In the definition of tlie possible breaking terms of the STIs and WTIs one first admits 
any kind of local Lorentz invariant terms which have the proper quantum numbers. 
However, the consistency conditions of Eqs. ( 3.12|) - (|3.14 ) constrain those breaking 



terms further. Actually, they play the key role in the algebraic analysis of anomalies. 
They single out the possible candidates for breaking terms which cannot be removed by 
suitable counterterms. It is well known that this can be done by means of cohomological 
methods. For this important issue we refer the reader to the rich literature ^3[|. As 



an example, in the SM the consistency conditions single out the Adler-Bardeen-Jackiw 
anomaly. However, as is well known, the latter cancels out because of the specific 
choice of the fermion content in the SM. 

Once we know that the identity has no anomaly, the practical constraints of the consis- 
tency conditions on the breaking terms (introduced by the non-invariant regularization) 
can be worked out. This is achieved by the condition that only those breaking terms 
are possible which can be produced by local counterterms. We will explicitly illustrate 
this practical procedure for the process 6 — 57 in Section |[ 

2. Another important consequence of the consistency conditions is the triangular structure 
of functional identities. This means that it is possible to organize the set of functional 
identities into a hierarchical structure in such a way that we can restore the identities 
one after the other without spoiling those which are already recovered — as explicitly 



explained in Section in the example of 6 ^ S7. 



We have to solve the following problems 



W[,){T'ct) = A'^(A), (3.15) 



where T'^t = Jd^xC'^^{x) with the conditions (CT)-(|3^)- 



We know that, in absence of anomalies, the breaking of the Wess-Zumino consistency 
conditions ( p.l4|) are solved by the counterterms 



W(\)(r'^%") = A'^(A) 

Therefore, by introducing those counterterms in the Feynman amplitudes, or equiva- 
lently in the action F': F' — F' — T^^^ , we have the new system 

^^'wir'-T'^/') = A'^(A)-W(\)(F;,^")=0, 
where the new breaking term A's is explicitly background gauge invariant because 

of ( prrq ) and ( prri 



W(\)(A'5) = W;,)(A^) - W;,)5^(Fe^") = W(',)(A'^) - S', (A'^(A)) = 0. 
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This means that, in order to restore the STI, we need only background gauge invariant 
counterterms. To that purpose, we add terms hke T^^^ so that 

W(',)(r;f^") = 0, S',{T'S^^) = A's. (3.16) 

This proves that we can effectively disentangle the two sets of identities. In Appendix |C| 
we specialize the argument to each subspace of counterterms. 

In the following we analyze the practical algebraic renormalization program within the 
BFM discussed in Section ^]3| in more detail. 

We consider the STI ( p.6|) and separate the extended part from the conventional one 
(namely the old STI (p.6|) ) by taking the derivative with respect to Q and then setting them 
to zero. As mentioned in Section |2.3| , we obtain a set of the following three different equations 
(in the following we regard only the example for but it is understood that this is valid 
for every Q): 

0, >V(A)(r) = 0, (3.17) 
T7^] (3.18) 




where F = F^^^q, iSr is the linear (conventional) STI operator ( |A.8| ) and W(a) is the WTI 
operator defined in (2.24). 

Notice that reducing the equation for an arbitrary f2 to = might hide some relevant 
Green functions with higher powers of Q. In particular in the Eqs. ( |3.17|) we consider only 
the first power of Q and the functionals independent of Q. However, we would like to recall 
that f2 is a vector or a scalar field from the Lorentz point of view, they have ghost number 
equal to +1 and dimension equal to 1. Therefore the higher powers of Q cannot contribute 
to divergent Green functions or, equivalently, there are no counterterms with two powers of 

n. 

At the quantum level they are broken by the breaking terms controlled by the QAP 

As + 0{hAs), W(A)(r) = A^(A) + 0{hAw) (3.19) 
- '^r (^] = Acs + OihAco^). (3.20) 




V^^M/n=o 

Notice, furthermore, that there is no prime in the following equations since we restrict our 
considerations to the functional F and to the breaking terms A5 and Aw in the case where 
f2 = 0. Finally we can use the consistency conditions (cf. Eqs. ( |3.12| )-( p.l4| )) and one can 
easily deduce the following set of constraints 

(3.21) 

(A5) = So (AgA (3.22) 



So {As] 


) = 




6 







29 



So (AwW) - W(A) (As) = (3.24) 

>V(A) (Agj) - - ^) (^^(^)) = fabcX'Aa^^. (3.25) 

This list is not complete because we have also to add the consistency conditions between the 
breaking terms Avi/(A), but these are exactly the same as written in (3.10). 

Using this equation we can study the two approaches singled out in Section (|2.3|) . In 
the first approach we restore the naive STIs, that means we compensate the breaking terms 
As with suitable counterterms. In addition we consider the WTI and their breaking terms 
A^(A) and we compensate them by using suitable counterterms. Those counterterms are re- 
lated by the consistency conditions ( p.24| ) which we take into account in our explicit examples 
in order to reduce the eff'ort of our computation. 

By setting all the breaking terms As = Avi/(A) = to zero and assuming that all the 
necessary counterterms are introduced, we have finally 



'^0 (Ag° 



(3.26) 



6 6 




^oi) - TAT - 77^ {^gA = (3.27) 



W(a)(Ag^) =/a6cA'AGc. (3.28) 



The first equation ( |3.26D tells us that the breaking terms A^c of the remaining identities 
(which are needed in order to control the breaking of the relations between the background 
gauge fields and the quantum ones) are STI invariant breaking terms. Therefore we need 
only counterterms which satisfy the STI. The second equation ( |3.27|) provides a consistency 
condition for the breaking terms A^c themselves. It tells us that the breaking terms can 
only depend upon the linear combinations of the background gauge fields. The last equation, 
namely ( p.28| ), says that the breaking terms A^-c should transform as vectors under transfor- 
mations of the background gauge symmetry. It is easy to see that the possible candidates for 
the breaking terms A^c are very few and they can be reabsorbed by adjusting counterterms 
with background fields. 



In the second approach we first restore all the WTI and all the identities like (|2.21 



That is we suppose that we have inserted the relevant counterterms which render A^c = 
Avi/(A) = 0. From the consistency conditions above one easily derives again the constraints 
on the breaking terms A5. 

In order to select a practical criterium to use one approach instead of the other we have 
to study the process under consideration. In fact two situations can in general occur. If we 
consider only processes at one- or two-loop, we have 

• the case where both, the quantum and the background version of a Green function, 
occur as sub-diagrams or as a component of the connected Green functions involved in 
the process under consideration. Here the second approach seems to be advantageous. 
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• the case where the quantum sub-diagrams are completely unrelated to the background 
Green functions involved in process (for instance in our example in Section \2.4\ ). Then 
identities (3.18) are not relevant and one is allowed to restrict the algebraic analysis 
to the WTI and to the conventional part of the STI (3.17). 



3.3 Removing the breaking terms 

In the following a general procedure is presented which optimizes the algebraic method. 



In ^J] this procedure was used to discuss the complete renormalization of the Abelian Higgs 
model. In this paper we use the procedure in the more general context of the SM. The 
breaking terms to be computed are reduced to the evaluation of finite Green functions. 
Moreover, it is shown that this method is very powerful if in addition the BFM is adopted. 

First a formal derivation of the various steps is presented. Afterwards the power of 
the method is demonstrated in the case of the two-point function for the W bosons which 
represents an important ingredient for the calculation of radiative corrections to the well- 



known p parameter |^ parameterizing the isospin breaking of a fermion doublet. 

Let us assume that the invariant vertex functions r*^™^ has already been constructed up 
to order m < n — 1. Thus, because of the QAP (see Section |3.1D , the subtracted functional 
r(") satisfies the broken STlQ 

n-l 

cS(r)(") = 5o(r(")) + ^ (r^-'), r^"-^)) = a^") , (3.29) 



where the meaning of Sq and the compact notation for (-, ■) is given in Eqs. (|A.8|) and ([A.9|) . 



In the case of the absence of anomalies, we know from the general results of algebraic 
renormalization that one can add finite counterterms in the action in order to restore the 
validity of the STIs. As we have seen in Section |2.1|, the STIs connect a large amount of 



Green function. In principle, all of them have to be calculated in order to fix the breaking 
term A^"). 

The construction of an efficient and convenient method for the determination of A*^") 
consists of the following steps: 

1. We disentangle the various terms in A*^") expanding it into an operator basis 

A(") = ^ Q y d'^xM^ix) (3.30) 

i 

where Aii{x) are local Lorentz invariant monomials in the fields and their derivatives. 
The basis is quite restricted by additional symmetries preserved by the regularization 
procedure. The usual power counting poses an upper bound on their mass dimension 
(which is independent of the loop order in the case of power counting renormalizable 
theories). As we have seen in Section A*^") is also constrained by the Wess-Zumino 



consistency condition given in ( |3.14|) 
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In the present section we skip the prime on the effective action again. 
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2. One acts with the zero momentum subtraction operator (1 — T^) on both sides of 
Eq. (|3.29| ). Here denotes the Taylor operator in the external momenta up to a 



suitable degree 6 (see Appendix ^ Eqs. ( A.IO ) for its explicit form). The locality of 



A*^") ensures that it is possible to find a degree 6n such that0 

(1 - T'^")5(r)(") = (1 - T^") A(") = , (3.31) 

and thus A^") is subtracted away. At the moment we assume that the zero momentum 
subtraction is possible. This means that the vertex functions have to be sufficiently 
regular at zero momenta. Below possible adjustments for the case that IR problems 
occur are discussed. 



3. Clearly the l.h.s. of ( p.31|) has not yet the correct form. Indeed we want to obtain a 



new STI for subtracted Green functions, i.e. for (1 — T^p'')T^'^\ where 5pc is the naive 
power counting degree. In general we have 5„ > 6pc. To that purpose we commute 
the Taylor operation with the Slavnov- Taylor operator. It is convenient to adopt the 
decomposition (|3.29| ) into a linearized operator plus bilinear terms. The part involving 



the linearized operator leads to 

(1 - T'^")5o(r(")) = So ((1 - T^f=)r(")) + So (r^r^cY^^ _ (5o(r("))) , (3.32) 

which expresses the fact that So is in general not homogeneous in the fields. In par- 
ticular this is the case for theories with spontaneous symmetry breaking. Notice that 
the Taylor operator is scale invariant, however, it does not commute with sponta- 
neous symmetry breaking. Furthermore there might be IR problems in connection 
with massless fields as will be discussed below. The difference between 5„ and 6pc leads 
to over-subtractions in r*^"'^ Therefore breaking terms generated by the last two terms 
on the r.h.s. of the above equation have to be introduced. Furthermore, the action of 
the Taylor operator can be split into the naive contribution E"=i (r^^), T^""^)) plus 
the local terms obtained by Taylor expansion. These local terms also contribute to the 
new breaking terms. 

Finally, by applying the Taylor operator on ( p.29|) and by using ( p.31| ) and ( |3.32|) we 
obtain 

n-1 

(1 - T^")5(r)(") = 5o ((1 - T^f-)r(")) + J2 {^^^\ r("-^)) (3.33) 

n-1 

- [T^^So - 5oT^*'=]r(") - T^" (^^■'^ ir^"^-''^) = . 

The terms in the second line of ( p.33|) represent the new local breaking terms which 
correspond to the subtracted function at the order n, (1 — T^p'')T^"'\ Thus, it is 
convenient to define the new breaking terms as 

n-1 

^(n) ^ ^T^r.^^ _ ^qT^p^JF^") + T^" (^^■'^ r^""-''^) . (3.34) 
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In power counting renormalizable theories the degree 6 is of course independent of the loop order n. 
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We emphasize that they are universal in the sense that they do not depend on the 
specific regularization used in the calculation — in contrast to A*^") in Eq. (|3.29|) . 

4. Now the principal construction of an invariant Green function F*^") is clear: First, 
one has to calculate the universal terms This step consists of the evaluation 
of a set of finite amplitudes and their derivatives at zero momenta as one can read 
off Eq. ( p.34| ). There, the functions r^™"-* are computed at the lower orders in the 
perturbative expansion. They are supposed to satisfy the STI at every order m < 
n. One strategy to further simplify the calculation of the universal breaking terms, 
^!'^'^\ is to reduce the number of bilinear contributions by a suitable choice of the 
renormalization conditions. Second, one has to find the counterterms S*^") which satisfy 

5'o(S(")) = . (3.35) 

Finally, the correct vertex function reads 

jp(n) _ _ 2-'5pc)r(") + . (3.36) 

We stress that only the universal breaking terms have to be computed which drastically 
simplifies the determination of the non-invariant counterterms. 

5. As mentioned above, in the presence of massless particles zero- momentum subtractions 
of the regularized function F*^") might lead to IR divergences. In principle one can 



circumvent this problem by using the Lowenstein scheme [0. In this case one has 
to introduce a generalized Taylor expansion which also takes care of the soft mass 
parameter s of the Lowenstein scheme. Consequently one has to analyze the new 
breaking terms arising from the commutator between this new operator and the STI, 
respectively, the WTI operator. Moreover, one can take advantage of the fact that the 
breaking terms are IR safe for principal reasons provided there are no IR anomalies in 
the model. In the phenomenological examples we discuss in the following sections no 
IR problems occur. 

There is also an alternative path for the case where IR problems are induced by zero- 
momentum subtractions. The whole procedure proposed in this section, that is the 
translation of the "conventional" breaking terms into "universal" ones, also leads to 
drastic simplifications if the subtractions are performed for infinite instead of zero 
momenta. The former procedure is even preferable if four- and five-point functions 
occur in the STIs. It obviously circumvents all IR problems mentioned above. 

6. We emphasize that there is the free choice of the renormalization conditions which 
corresponds to a specific choice of the invariant counterterms. A change in the renor- 
malization conditions leads to a change of the basis of the breaking terms. In practice, 
this means that one is allowed to shift the breaking terms of some WTIs and some 
STIs to others which are not relevant for the specific process under the consideration. 
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7. There are also cases where the zero momentum subtraction is not very practical. Let 
us consider a STI (or a WTI) involving Green functions with four external legs as a 
minimum (e.g. Tw^w^ApA^) or with high- dimensional fields like the BRST sources (e.g. 
^c+q*b where q* is the BRST source for the quark field g,). As is easy to show they are 



fixed by STIs involving convergent Green functions (see Sections |2]^ and [5.2| ). There- 
fore zero-momentum subtraction introduces breaking terms which are very tedious to 
compute as they in general involve Green functions with five or six external legs. Very 
often quite a lot of diagrams contribute. Thus sometimes it is more convenient to 
choose another subtraction technique. One possibility would be to consider large ex- 



ternal momenta. Using Weinberg's theorem |3^, the coefficients of the breaking terms 



are computed directly by Green functions with only a small number of external legs, 
Uke Tc+q*b- 

8. The subtraction technique for STIs and WTIs correlated to Green functions like Tc+q*b 
or similar ones are more involved as can be seen from Eq. ( |2.55| ). Thus it is convenient 
to decompose the STI ( |3.29|) into different terms: 



n— 1 n— 1 

j=l a^U a£U' j=l 

= A("), (3.37) 

where U, respectively, L(' are index sets corresponding to the n-loop Green functions 
r^") and r^"). Remember that the finite parts of F^") still have to be fixed whereas this 
is already done for r^"\ The functions F^") are conveniently computed using other 
techniques, e.g. the one described in 7. The coefficients Cq and da are functions of the 
external momenta. 

Now the Taylor operator (1 — T^) can again be applied leading to the following result 

n-1 

^ ca ((1 - r'^-)ri")) + ^"iri") + ^ (r(^), r("-^)) + (3.38) 

I3eu aeW j=i 



n-1 



aeW 3=1 







where the new contribution T^" ^am' ^a^^a^ appear in the breaking terms \1>("). It 
depends on the Green functions F^"-'. This modification of our formalism will be used 
in Section ^ for the example h 57. 

In order to illustrate the method elaborated above we consider the two-point functions 
for the W boson. All the ingredients to consider the STIs involved in the calculation of 
the two-point functions in a non-symmetric regularization scheme are now introduced. In 
Appendix ^it is shown that the two-point functions and their STIs (plus the ghost equation) 
form a closed set of relations. It is also known that the STIs (|2.48|) are broken by local terms 
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only. According to the QAP we have: 

^c+w;''iP)^w+w-iP) + ^c+G*^-ip)^G+w-iP) 

W+G-iP) + ^c+G'-''{p)^G+G- ip) 



A,+G-{p)+0{nA). (3.39) 



The breaking terms A^+i^^-(p) and Ac+g-{p) are finite and exhibit the following decomposi- 
tion in terms of independent monomials 



A 



w-iP) 

^c+G-iP) 



aiP Pfi + a2P^i 
asp"^ + . 



(3.40) 



According to our procedure we apply the Taylor operator to the various Green functions 
involved in the STIs 





(p) = 


(1 


rTl2 

p 


)^W+W^iP), 




^G+W- 


(p) = 


(1 


p 


)^G+W-iP) = 


-ip^(l - T^)Tg+w-{p'') , 


^W+G- 


(p) = 


(1 


rp2 

P 


)^W+G'iP) = 


-ip^{l - T^)Tw+G-{p'') , 


Tg+g- 


(p) = 


(1 


rp2 

p 


)rG+G-(p) , 




f 


(p) = 


(1 


-T' 

p 


)^c+W*'~iP) = 


-zP^(1-T;)J(p^) = -tp,J{p^ 


Tc+G*.- 


ip) = 


(1 


-T' 

p 


)Te+G'-iP) = 


(i-T,W) = V), 



with r 



(3.41) 

c+w*~(p) ~ ~WvJ{p'^) and Tc+g*-{p) = lip^)- The UV index is set equal to the 
superficial UV degree. 

In a next step we can compute the commutator between the Slavnov- Taylor operator 
and the Taylor subtraction. This leads us to the universal breaking terms (cf. Eqs. ( |3.33| ) 
and dpI)) 



(1) 

c+W- 



f^>(i) ^J2^ 
w+w- 



\P' 



iMwT^gIw- 



{p')-zM'^j('\p')-MwP\p' 



tp^p 



d 



-^Mi^— rg|^^_(0) - Mh./(^)'(0) - ^M^J«'(0) 



(3.42) 



(1) 

c+G- 



P '- W+G- 



{p') + tMwr'ilG-ip')-p'MwJ^'\p')+^p'i^'\p' 



0. 



where AT^q+w- 



(0) 



dp^^ G+W' 



ip' 



p=0 



w+w- 



denotes the longitudinal part of the W 



boson two-point function defined through 



^w+w~iP) 



p^ 



{P' 



In the second equation of (p.42|) there is no breaking term at all while in the first one 
only the term proportional to p^p^ survives the Taylor operation. It can be re-absorbed by 
introducing the counterterm 



MH.(-^rg|^40)+2/«'(0)-M^J«'(0)) j d'xd^W;:d''W; , (3.43) 



dp 
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which is manifestly universal because it is computed in terms of finite Green functions which 
do not require any regularization procedure. This is the major difference between the direct 
computation of the breaking terms and the one discussed here. In the former case the 
complete calculation of the coefficients ai, 02, 03 and have to be performed. On the 
contrary, in our approach this is in part automatically taken into account by the subtraction 
procedure. In the present example, we only have to compute the residual coefficient ai which 
arises in Green functions involving ghost fields and W-G mixing. The other breaking terms 
of Eq. (|3.40| ) are automatically removed by the specific choice of renormalization conditions 
given by the Taylor subtraction in (|3.41|) . In other words it is possible to remove the breaking 
terms {i = 2,3,4) by invariant counterterms. 

Note that in Eq. (|3.43| ) only the first term in the bracket gives a contribution to the 
fermionic part since at one-loop order the ghost two-point functions do not receive any 
contribution from fermionic fields. 

Before we discuss the additional simplifications of our method in connection with the 
BFM, we mention that our specific procedure allows us to work out the corresponding 
formulae at the n-loop level explicitly: 



c+W, 



n-l •] 

^^_(0)-Mh./(")'(0)-^M^JW'(0)] 



m=l 



^P.{p' i-^^^^^cl 



n-l 



-p'Y. 



m=l 



j(-)'(o)r^(r^-)(o) + :r(")(o) Ar^(— )(o) 



+ j(-)'(o)r(,v™l (0) +x("^)(o)^rJ,-;^) 



(0) 



(n) 



c+G- 



n-l \ 
m=l ) 

n-l 

+ E + j('")(/)r(,VcJ(p) 

m=l 
n-l 

= E [-p'j^"^^ mt^t (0) + x(™) m^G^G^ (0) 

m=l 



m) 
G+G' 



(0) 



(3.44) 



where -ip^j''''\p'^) = r|_^^^.,_(p) and l^^\p^) = r^+g„,_(p) (A; = - 1) represent 

the renormalized symmetric counterparts of r|^^^^,,„(p) = —ip^jj''^\p^) and Tf}^^^^{p) = 



36 



l{k)(jP'y Already at two-loop order the complete set of counterterms is needed. They are, 
however, expressed in terms of one-loop Green functions and only for Oi a real two-loop 
calculation is needed. The closed formulae for the coefficients of the breaking terms, ai'\ 
^2"^ c^s"^ and 0:4"^ where the superscript "(n)" reminds on the number of loops, allow us to 
adopt special renormalization conditions which simplify the expressions. The choice 

J^^\id) = l, X(°)(0) = Mh., :r('=)(0) = 0, J('^)(0)=0, k>l, (3.45) 

in combination with the adjustment of the wave function renormalization of the ghost and 
Goldstone field, Zf and Z^, reduces the breaking terms to 



= ^P.{/(-^^^^r[,1^_(0)-MH^/(")'(0)-zAf^j(")'(0)' 



n-1 



-p^;^ [:r(-)'(o)r^?,^-)(o)+j(-)'(o)rg^-™)(o)^ , 

m=l ) 

n-1 



m=l 

„W ov.^ ^,(") 



where the coefficient 0:2 and 04 are zero to all orders. Therefore, only the counterterms 



J d^x [at'^d''W;ld''W- + af^d^G+d^G-) , (3.47) 



are needed to restore the symmetry. The appearance of the lower order terms is a spe- 
cial feature of the STL In fact this is strictly related to the renormalization of the BRST 
transformations at the quantum level. 

Finally, we show that the BFM in combination with our subtraction procedure provides a 
further simplification in the S-matrix calculation. Instead of using the conventional approach, 
the S-matrix elements can also be computed by the BFM. In that case the above identities 
are replaced by 

¥^l^. = p2f«^G-(p)+^^w^r«^_(p) = 0, (3.48) 

where and ^^^^c- are the universal breaking terms of the WTIs ( |2.25| ) and ( |2.26| ), 

respectively. We stress that the same WTIs hold to all orders without any changes. This 
is due to the linearity of the identities. No ghost fields occur in these equations. As a 
consequence, only the counterterm 

S(^) = (-^r«^_(0)) / d'xd^W;:d'^W- , (3.49) 

has to be introduced in order to restore the background gauge symmetry. We mention that 
the corresponding counterterm at the order n is given by the same formula. 
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At this point we again want to underline the important difference between the BFM 
approach and the conventional one. In the former case, the computation of the amplitude 
can be performed by using Green functions with external background fields instead of gauge 
bosons and scalars. For these Green functions only WTIs are needed to define their finite 
parts correctly. However, at higher orders also quantum two-point functions appear as sub- 
diagrams. This implies that STIs have to be employed at lower orders. By using the method 
exposed here we can immediately see that for one-loop STIs already the zero momentum sub- 
traction partially restores the identities while only few additional counterterms are needed. 
Furthermore the restoration of the WTIs is simple anyway. The method of subtractions 
provides a very powerful technique to renormalize a model in a non-invariant regularization 
scheme, especially if in addition the BFM is adopted. 

In the forthcoming sections, we will see that in the case of an external photon field (as 
in 6 — i> S7 and H 77) there are further simplifications due to U{1) group of the SM. 

4 Two-loop QCD corrections to H ^ jj 

In the minimal version of the SM the only not yet discovered particle is the Higgs boson. 
Up to now only limits on its mass exist where the lower one is provided by the direct search 
at LEP. The upper limit is obtained from the combination of precision measurements and 
quantum corrections. A Higgs boson in the intermediate mass range, i.e. Mw < Mh < 2Mw, 
is very promising. Then the loop-induced decay into two photons is one possibility for its 
detection. Next to the diagram involving gauge bosons also the one where the Higgs boson 
couples to top quarks accompanied with additional QCD corrections is important. In the 
latter case an expansion for a heavy top quark mass is sufficient in order to obtain a sensible 
approximations to the exact result. In the following we will exemplify our approach to the 
algebraic renormalization developed in the previous sections in this limit. The calculations 
are performed in the framework of Analytical Regularization which will be introduced in the 
first subsection. 



4.1 Analytical Regularization 

Although quite a lot of different regularization prescriptions have been developed since the 
invention of quantum field theory by far most practical applications have been performed 
in the framework of Dimensional Regularization p, |^ . Meanwhile quite some technology is 
available which allows the computations of rather complicated Feynman diagrams sometimes 
even at the four-loop level. One of the main reasons for this development is that most of the 
symmetries are preserved. Furthermore, very often the occurring d dimensional integrals are 
even simpler to evaluate than the four dimensional ones. 

A regularization method which was invented even before the advent of the dimensional 



one is called Analytical Regularization (for the practical details we refer to |2^). In the 
dimensional method the space-time over which the integration is performed is changed from 
four to a complex number, d, and the limit d — >■ 4 is taken at the very end. The divergences 
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then appear as poles in l/{d — 4). On the contrary in the case of Analytical Regularization 
the space-time is kept fixed and only the exponents of the denominators are modified. One 
of the big advantages of this approach is that no problems in connection with 75 (or, more 
generally, the e tensor) appears. The integrals which have to be solved can be slightly more 
difficult. It is, e.g., not possible to perform a partial fractioning as the exponents are no 
integer numbers any longer. For some classes of diagrams it is nevertheless possible to take 
over the results obtained in Dimensional Regularization. As an example we consider the 
massless scalar one-loop integral. For d 7^ 4 and arbitrary exponents the result reads: 

d'^q 1 
{2nY {{q -pY + ir])''{q'^ + ir]y ~ 
(_p2)d/2-a-6 r{a + b-d/2)r{d/2-a)r{d/2-b) 



(4vr)<^/2(_i)-a-6 T{d-a-b)T{a)T{b) 

This result can immediately be interpreted in the framework of Analytical Regularization 
after choosing d = 4 and a = ai -\- X and b = bi + fi where ai and 61 are integers and A 
and /i adopt the role of the regulator. Possible divergences appear in the form 1/A, and 
1/(A + /i) after the r.h.s. of Eq. ( [4.1| ) is expanded in a Laurent series. This also leads to 
logarithms with dimensional arguments. Therefore, one has to introduce an arbitrary mass 
scale, fi, such that this is corrected for. 

Dealing with different regulators for different denominators would be in general quite 
tedious. In our case the structure of the divergences of the diagrams allows us to choose 
non-integer exponents only for the gluon line — the fermion lines have not to be regulated 



at all |4|. 



4.2 Breaking terms for i7 ^ 77 



According to Section p.2| we have to discuss the structure of the breaking terms of the STIs 
for the amplitude ( p.l8| ) and the sub-divergences ( 2.15| ). Recall that we are only interested 



in the two- loop QCD corrections to the fermionic contributions IT^^^^^^. 

The most general Lorentz invariant breaking terms corresponding to ( |2.15| ) and ( |2.16| ) 
are given by 

Ac^aM^^^p) = (5!'^ + 5f0^2"+(5^'^ + 5r)p^ (4.2) 

where and are the one- and two-loop coefficients needed for the vertex contribu- 
tion IT^^^ {i = 1,2). 6^^^ and ^4^^ are the sub- leading breaking terms for the sub-graph 

The breaking terms are removed by expanding the Green function T^HApA„ {lii ^2) w.r.t. its 
external momenta and discarding the first non trivial contribution as the operators {1—T^_^ 
have to be applied. Here T^^^^^ denotes the Taylor operator up to first order in qi and q2- 
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For dimensional reasons the constant term ^4^^ has to be zero. The computation of 
Sl^^ is discussed in the forthcoming subsection. Again it is enough to apply the operators 
(1 — Tpp) and (1 — T^) to the corresponding Green functions, r^^'^^*^^ (p, p) and T^^l'l^'^^ (p), 
respectively. This corresponds to adjust the free parameter Zi in Eq. ( p. 131) . 

4.3 Results 

Since we restrict ourselves to two-loop QCD corrections, only diagrams involving the top 
quark have to be taken into account. Some sample diagrams are shown in Fig. ^ For our 
purpose it is enough to consider the limit where the top quark mass is much larger than 
the Higgs boson mass. This allows us to perform an expansion of the vertex diagrams in 
the external momenta and thus reduce the integrals to be solved to vacuum graphs. We 
note that the result is known in Dimensional Regularization [BB, ETl which allows us for a 
comparison at the end. 

It is convenient to split the contributions to the H'-f-f vertex according to the loop order 
in the following way 

^HAf^A-iqi, 92) = r^AMA-l^l. (I2) + — CFr^AMA-l^l. ^2) + • • • , (4.3) 

TT 

where the momentum and polarization vector of the photons is given by (gi, e'^) and (g2, e*^), 
respectively. The general Lorentz decomposition of the vertex diagrams is given byQ: 

Note that in regularization schemes where gauge invariance is preserved — like Dimensional 
Regularization in the absence of chiral fermions — one has according to ( p^.lSj ) A = —B. 



In the case of Analytical Regularization this is not true and some breaking terms occur. 
However, it turns out that these breaking terms can be removed by Taylor subtraction. 
No new universal breaking terms have to be introduced as in the case of over-subtraction. 
Therefore the subtracted Green function (using Analytic Regularization) fulfills Eq. (|2.18| ) 
at the one- and two-loop level 

-zg^'(l-T° Jrg^^^^(gi,g2) = 0, 2 = 1,2 (4.5) 

where qi and q2 are the momenta of the photons. This relation can also be rewritten in 
terms of the functions A, B and C with accordingly adjusted Taylor operators. 

Let us now consider the one- loop calculation, which constitutes the Born result, in more 
detail. The evaluation of first three terms in the expansion for large top quark mass leads 
to: 

At 

= A(— + l + —f+—f^ + Oif')], (4.6) 
V4r 30 21 ^ V ^ ^ 



14 



In the limit of on-shell photons there is no contribution from C to the decay rate. 
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with A = NcQf2a/ {3ttv) and r = 2gig2/(4m^). rrit is the top quark mass in the MS scheme. 
As Eq. (|4.6| ) constitutes the lowest order, may as well be replaced by the on-shell mass, 
Mt. Clearly the application of the operator (1 — T^-^^g^) to qi.q2A^^^ removes the breaking 

term proportional to 1/f and leads to the invariant Green function F^^^^^. 

At two-loop order the situation is different. Here, the divergent sub-diagrams have to be 
carefully investigated in a first step in order to fix the renormalization constants Zi, Z2, 
and Zy defined in Eq. ( |2.13| ). As described in Section ^.2\ , it is possible to fix Z2 and Z^ in 
the MS scheme and to choose Zy = Zm- From the examination of the fermion propagator 
we get 

TT 8A 

Z„ = Zy = l + ^C, (-^) . (4.7) 

Here, ^ is the gauge parameter appearing in the gluon propagator defined in the conventional 
way i{—g^'^ + ^p^p'^ / p^) / {p^ + irj). Thus we are essentially left with Zi. Its divergent and 
also finite part is fixed by the QED- WTI given in Eq. ( p.l5| ) . We again use Taylor expansion 
to get 

^{P + pT (1 - Tl) T^ZfriP^p) + -Q. (1 - [r^t^'^""(p) - ra^'«""(-p)] 

= 0, (4.8) 

where p and p in the three-point function denote the momenta of the quark and anti-quark, 
respectively. Note that also here no over-subtraction is introduced. The Taylor coefficient 
'^P,p^Al^^q^^iP^P) fi^^s -^1- turns out that the finite part of Zi is zero and the divergent part 
fulfills the WTI. This means that the Green functions obtained in Analytic Regularization 
fulfill the WTI even without Taylor subtraction. Thus we have 

^1 = ^ + -Cf^^. (4.9) 

TT OA 

At this point it is convenient to introduce the notation: 

gi.g^A^'* = (l-T°,^Jgi.g2A», (4.10) 

and in analogy for the functions B and C. Then the renormalized function A up to 0{aas) 
can be written in the form 

A(q,,q2) = (1 + + 25Zi - 35Z2) + — 



+ ^^^(2),*, (4.11) 
m'^=Zmmt 



where m° is the bare top mass and the condition Zy = Z^ has been used. In the first line 
the first occurrence of bZ2 origins from the Higgs-quark vertex. Furthermore, all three quark 
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propagators deliver a 6Z2 and each photon-quark vertex leads to a 6Zi. After the second 
equality sign SZi = 6Z2 has been used. 

In Analytical Regularization the two-loop result reads in the MS scheme: 



A 



3 
1 



181 



360 20 



/1541 



mr 



V 4725 7 



f 2 + 0{f 



mi 



. (4.12) 



In order to compare with the results known from Dimensional Regularization physical 
conditions have to be imposed. We choose to express the result in terms of the on-shell mass, 
Mt. It is defined through the zero of the inverse propagator where the external momentum 
is on the mass shell. In Analytical Regularization the transformation from the MS mass to 
the on-shell one reads: 



mt 




3 3 , /i 
m 



The substitution of this equation in (|4.6|) leads to 

3 61 2216 



A 



4 



Finally, the decay rate is given by 

T{H ^ 77) 



90^ ^ 4725 



4 "' Mf 



(4.13) 



TT 



Ml 



H 



647r 



(4.14) 



(4.15) 



where the functions A*^*) have to be evaluated for qi.q2 = M'jj/2. The decay rate T(H 77) 
coincides with the result obtained in Dimensional Regularization |^6|, 

Although the result presented in this section is quite simple, the main steps of the Alge- 
braic Renormalization have been touched. We have also seen that in this specific example 
our algebraic procedure is as efficient as Dimensional Regularization due to the linearity of 
the STI involved. 



5 Two-loop electroweak corrections to 6 ^ 57 

In this section, the electroweak two-loop corrections to b ^ s'j are discussed. The inclusive 
mode of this rare decay is already measured and plays an important role in the search for 
physics beyond the SM (see e.g. ||38|). In the literature a partial calculation of the two- loop 
electroweak corrections has been performed |3^. In the limit of a heavy top quark and/or a 



heavy Higgs boson also a complete calculation exists |]31[. In all cases the naive dimensional 
scheme for 7^ has been used. 

In our analysis, we use the BFM in combination with our practical algebraic framework. 
All relevant WTIs and STIs are explicitly derived in Section [^.4|. Special care is taken in 



order to explicitly illustrate every step of our subtraction method presented in Section ^ 
A comparison with the conventional algebraic method is made. 
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In the following A^'^") and A'^''^") denote the n-loop breaking terms to the WTI and 
the STI, respectively. For the differentiations w.r.t. fields it is convenient to introduce the 
notation: 



^n^W/5,(n) 



6(f)l{pi) . ..6(f)niPn 



^W/5,(n) 



0=0 



1) 



(5.1) 



5.1 Restoring the WTIs 

Let us in a first step consider the WTI of Eq. ( p.33|) 

5^>V(A)(r)(2) 



5XAi-Ps - Pb)Ss{ps)6b{pb 



KPs+PbT'rf^^{Ps,Pb) + ieQa 



.(2) 
Sb 



<f>=0 

ipb) 



sb ^ 



-Ps) 



^\Asb 



iPs,Pb) 



(5.2) 



The corresponding counterterms (p.34|) have already been presented in Section |2.4| . They 
are fixed by hermiticityf^ and Lorentz invariance. Furthermore they are constrained by 



the consistency conditions given in the Section |3.2| . In order to impose them we follow 
the suggestion presented there. Since we know that in the SM all breaking terms can be 
removed by non- invariant counterterms or, equivalently, no anomalies are present, we insert 
our general ansatz of Eq. ( p.34|) into ( |5.2|) and deduce the most general breaking terms 
compatible with the consistency conditions. At the two-loop level the most general form 
reads: 



W,(2) 
^XAsb 



{Ps,Pb) 



Pl 



4 + alfPR 



'Is + a'sb^Pn 



(5.3) 



Actually the coefficients 13].^^^ are zero and a'^l'^^^ = —a]^^^^ as can be seen by using the 



1,L 



consistency conditions in combination with the invariance under charge conjugation. a^{^ 



and a^f,^ are fixed by (|5] 



a 



l,L/R 



R/L 



ip .Jd ^^W(.)(r)(^) 

4 ^/^^ \dp-^ 5\A{-Ps-Pb)5-s{ps)6h{p,) 



iPr/lI 



rJ>,0)-eQ.- 



d 



sb \ Ps 



<f>=0. 



P8=0J 



Ps=Pi, = 



(5.4) 



This equation shows that at the two-loop level the corresponding finite counterterms are 
identical to the first terms of the Taylor expansion of the Green function. Therefore, by 
using our subtraction method we automatically restore — up to sub-divergences — the WTI 
at two loops. 



One can eventually choose non-hcrmitian currents. Then, however, the number of counterterms in- 
creases M. 
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Let us move to the discussion of sub- divergences. In the conventional algebraic method 
one starts again with the analysis of the breaking terms to the various WTIs ( p.33| ), ( p.38| )- 
( p.41| ) in the same manner as above. For example considering Eq. ( |2.39|) one gets the following 
general breaking term: 

+ (^6 QfiuP"" -p"" + 5^ g^yp' -p' + 5s g^,u + e^^p^p'^'^p''" , (5.5) 

where are the momenta of the bosons. From hermiticity one obtains ^3 = ^4 and 
5q = 5j. Using the consistency conditions, i.e. inserting ( |2.47D into Eq. ( ^.39| ), leads to 
5^ = 0, + 25e = and 81 + 82 + 28^ = 0. Therefore, the r.h.s. of Eq. ( p.5|) reduces to the 
following four independent structures 



^'xaW+W,' (P^^P ) = "2 {p^pt - 2p+P^ +Pf,Pu)+(^3 {p^Pt - '^Pf.Pt + Pf^Pv 

+ ai9^.v - 2p+p-''' + p-pP''") + a^e^^p^p+'Pp-'"" . (5.6) 

The other breaking terms could be treated in analogy leading to the results 



^^agIw-^p'^^p') = 

^7aG+G- iP^^P~) = "8 {P^P^''' - 2p+p~''' + p~p~ 

= <.,PL(^.-^,) + «fo,,P«(^.-^,) , (5.7) 

where a^j^, ■ ■ ■ , C(]j'^ are matrices in the flavour space for u-, respectively, rf-type quarks. 

Clearly, the calculation of all the coefficients in Eqs. ( |5.6| ) and ( |5.7| ) involves all Green 
functions of the various WTIs and is quite tedious. However, a closer look to these equations 
shows that already a simple Taylor expansion of the Green functions can remove all possible 
breaking terms. Therefore, the Taylor subtracted Green functions (up to the order of the 
superficial degree) 

^iJp^^Ps) = i^-Tljrf^,M,p,), 

^'Zw-w.iP'-^P') = (^-Tkp~)^'Zw-w,(P^^P')^ (5-8) 

and analogously for the functions tf^^^^_{p+ ,p~), tf^^^^_{p+,p-), tf^^^+^^{p+,p-), 
^aIc+w;--^P~^'P~^' ^bJiPb,Ps), i'u^u.iPhPj) and t^^^^^{pi,pj) automatically satisfy the WTIs. 
Thus, they coincide with the symmetric Green functions denoted by T^^^ ^_{ph,ps), . . 

^TdM^Pj)- 
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At this point we mention again the subtleties involved in this procedure: if the degree of 
the Taylor expansion (necessary to remove the breaking terms) is larger than the superficial 
divergence of the corresponding Green function, local over-subtraction will produce new 
breaking terms to the WTI. Due to the fact that the WTIs (^) , (|]38D-(|^ do not 
explicitly depend on mass terms, no over-subtraction is introduced in this specific case. 
However, we will encounter this problem in the next subsection. 

The final result of the above analysis is that our prescription in combination with the 
BFM is a very powerful technique to get rid of the breaking terms. In particular no breaking 
terms for the WTIs must be computed and only an overall subtraction has to be performed. 
Regarding the WTIs with an external photon our procedure is as efficient as an invariant 
scheme since in this case simple identities (cf. Eqs. ( ^.2|) and (|4.2|) ), like in QED, are obtained. 
However, we stress that the finite counterterms corresponding to the first terms of the Taylor 
expansion have to be taken into account at higher orders, i.e. in our case the finite one-loop 
counterterms enter the two-loop calculation. 



5.2 Restoring the STIs 



Note that the STIs for the two-point functions (cf. Eq. ( p.48|) ) have already been discussed 
in Section |3.3| within our technique (see below Eq. ( p.41| )). The universal counterterms in 
the latter approach are given in Eq. ( p.43| ). Thus, in this section we only have to consider 
the remaining identity of Eqs. ( 2.52|) . 



In the conventional algebraic method, one would analyze the most general breaking terms, 
^'c+^^biPg^Pb) ^t-sqSPsyPq)' of these STI (|2.52| ) and the corresponding one obtained by 
the obvious replacements. Then one has to demonstrate how they are restricted on the basis 
of the consistency conditions by inserting the most general counterterms into the STIs. One 
could calculate these breaking terms in terms of the Green functions appearing in ( |2.52| ) by 
projecting out the different momentum structures of the breaking terms. 

However, the latter information is not needed if again our Taylor subtraction method 
is applied. 



S (1) 

Let us in the following explicitly consider A^iy^(pq,pfc). Similar equations hold 
for A^lg' {ps,pq). The family index is suppressed whenever there is no source of misunder- 



standing. According to the formula (H] 
respectively to ( p.53|) 



we apply the Taylor operator (1 



' Pb,Pq' 



to (P3^ ) 



^1 ^ 

Pb^Pq' 



(1) 



0=0/ 



Pb,Pq 



)A%)(p„p,) = 0, (5.9) 



where the degree of the Taylor operator is chosen to be the lowest one which cancels the 



breaking terms as is explained in Section 3^. Let us in the following explicitly write down 
the individual terms and subsequently discuss how they are treated. Eq. ( |5.9|) in combination 
with (|233| ) leads to 

= (1 - Tp\ pjr^J.^^.,_ + pb) (iVgbluPL) +i{Pq+ Pb)^ (1 - Tp° ,pjrJ^Vg-fe(Pg, Pb) 



+ - T^b,Pq)^%.APg + Pb) 



qb 



w 
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(P<7, Pb) 

A - m,)T%,,ip„p,)] - (1 - )r«(-p.) HV,,,Pl) 



V,,,Pn) (1-T^\)r«(p, 



Pb,Pq' 



^^c+qb* {Pq,PbWb-mb) 



(5.10) 



where in the second hne we have exploited that the first derivative of the scalar function 
rJ:+^,,_ vanishes. In the second term of the first line the operator T^^^^ has already been 



commuted with the momenta. Eq. ( 5.10 ) tells us that in most cases the superficial UV degree 
coincides with the one necessary to remove the breaking terms. Let us in the following con- 
sider those terms in more detail for which the Taylor expansion provides an over-subtraction. 



^G+qbiPi^Pb) has a logarithmic superficial divergence and thus we have 



G+qb 



(1 - T^.pXGiqbiP.^P'') = (1 - Tljr^GiqbiP.^Po) + K,pq - TU^^'ii.M^P^) , (5.11) 

where the last term corresponds to an over-subtraction. Also the terms proportional to the 
equation of motion for fermions lead to over-subtractions 



rpl ^ 

Pb,Pq' 



.(1) 

c+qb 



1 - ^p\,pJrl+V(P9>Pb)- - Tl,p)T^^,i,*iPq,Pb) 



(5.12) 



as T^^+-fj,{pq,pb) is only logarithmically divergent. These over-subtractions lead to the new 



universal (regularization independent) breaking terms 

]{Pq,Pb) = 



5,(1), 
c+qb ' 



^5,(1) 
c~sq 



{Ps,Pq) 



(TLpq-^Lpq) {MwT^GlqM,Pt 

+ ^q^^c!^rb(Pl' Pb) + ^bT^^gb' {Pq, Pb) 
{^Lpq-TLpq) {MwT^GKqiPs.Pq) 



+ msT^^^s„*{Ps,Pq) + rrigT^^}- {ps,Pg 



(5.13) 



where the Green functions r^+-*., r^+- ,. , F^j*. . and rj:j'-*^. will be discussed below. Inspec- 
tion of the most general counterterms available in the STI ( ^.52| ), namely those provided 
in Eq. (|23l| ) 



^(1) 

^W+q,b 
^G+q^b 
""qib 



Zi^'iPh + Zi^^Pr , 

Zi,5PL i>b + ^i,6^_R i>b + Zi^jP^ + Zi^^PR , 
Zij^dPl '^q + Zij^ioP^ + Zij^iiPi + Zij i2PR 



(5.14) 



(also r^^^ 



c+W, 



. _ and r^+^* _ discussed in Section R| should be taken into account) shows that 



S (1) 

the universal breaking terms, \E'^|y^(pg,pfe), can be removed by adjusting the universal coun- 
terterms and (compare with Eq. ( |3.35| )). Analogously '^'^-^^,{ps,Pg) can be re- 



moved by E^^__^ and Note that the Green functions ^c+g-biPqiPb) and ^G-sg.{Ps,Pq) 
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enter the two-loop calculation of 6 — 57 as sub-diagrams which means that they have to be 
computed anyway. Therefore, in this respect it is no extra effort to compute the relevant 
counterterms of our subtraction scheme as compared to an invariant one. 

However, in the calculation of the universal counterterms 3^^+.^, S^^l-^,, and 



-r^-o^ also the Green functions Tc+o*b-, 



\+q,b*, ^c-sq- and Tc-s*q,, appearing in ( p.l3|) . 



are 

needed up to linear order in the external momenta. In Section |2.4] , the complete set of 
identities for b —>■ sj has been discussed and Eqs. ( |2.54| ) and p.55| ) have been derived which 
constrain these Green functions occurring in our STIs. The WTIs ( 2.54 ) and ( p.55| ) are 
possibly broken by local terms like A^'^+-.^. 

One finds that within any mass-independent regularization scheme the breaking terms 
^r'c+o'fe ^'^'^ ^A^c+o*fe ^^^^ ^^"^ additional adjustment is needed. In the following, 

A Hi Z Hi 

we will explain the latter point in more detail. In addition, we also cover the case where a 
general regularization scheme is used. 



The breaking terms A 



and A 



Xzc+q*b 



could again be removed by Taylor subtrac- 



tion, however, over-subtraction would generate the following universal breaking terms (cf. 
Eqs. (1^ ) and (g3§)) 



X^C+q'bi-PQ-Pb,Pq,Pb) 



^\zc+q*b^ 



-Pq-Pb,Pq,Pb) 



0, 

-^^rgoc+,76(o,o,o) 



(5.15) 



where Mz is the Z boson mass and Gq is the corresponding would-be Goldstone boson. 



Note that ^clc+q'bi^^^^^) finite and could in principle be replaced by r^^^+-,^(0, 0, 0). 
For brevity the corresponding equations involving A+ have been omitted. 

As discussed in Section ^l3| , we will use our modified subtraction technique to compute the 
finite part of the remaining Green functions. In conjunction with the prescription provided 
there we found a very convenient and definitely more effective way to compute the breaking 
terms Ay^'^+-,, and Ay^'^+-., where all external momenta are considered as large: 

A Hi Z Hi 



(1) 



A 



W,{1) 



lim 



\AC+q;h 5Aa(0)5c+( 



aW,(1) 
^\zc+q*b 



-pPq - PPb)Sq*{pPq)6b{pPb) 

5^w{i;(r) 



ins. 6XziO)6c+{-ppq - ppb)6q*{ppq)5b{ppb) 



(/)=0 



(t>=0 



(5.16) 



For high values of the momenta the Green functions like F 



(1) 

Goc+q*b 



in ( |2.55| ) tend to zero 



because of Weinberg's theorem |32 



Also in the other WTIs we find that Ay^'*-!-*-*. and 

Aac+i?; b 



^A^c+o*6 combinations of the Green functions like F^+g*;, and T^+q.^ 

Z Hi i 



Note that the 

latter are simple to compute at one-loop order. Furthermore only one diagram is involved. 
Thus, the calculation of Tc+q*b and Tc+q^b* at zero momentum which is needed for Eq. ( |5.13| ) 
requires this simple additional step. Using this procedure Eq. (|3.37|) can be applied in 
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order to find the breaking terms \l/^+y^(pq,pfe) and Pg) where the symmetric Green 

functions ^c+q*b and r^+g^fc* can be identified with the "d-terms" of (|3.37|) . 

Note that in a mass-independent schemes, hke for instance Analytic Regularization and 
Dimensional Regularization with 't Hooft-Veltman 75, one shows that these breaking terms 
are zero, as already stated above. 

Unphysical renormalization constants as the wave function renormalization and the renor- 
malization of gauge parameters can be used to optimize further the algebraic technique as 
explained in Section ^. 



5.3 Imposing renormalization conditions 

As a last step, we have to fix the physical renormalization conditions in order to obtain 
the correct amplitude for 6 — ^ 57. We stress once more that the renormalization conditions 
must be imposed on the renormalized Green functions F by using the invariant counterterms 

4'i'sT[0 0[Il]- 

For the W boson we choose the following condition]^ 

rWf„-(p)]„.„,,.„. =0 (5.17) 



on the quantum self energies although it might not be convenient to fix the mass counterterms 

w+w- 



by using the conventional quantum two-point functions r^^V,!/- within the BFM. There, in 



fact, the renormalization condition on background two-point functions r|,\+,},_ are used. 

However, in the decay b ^ S'j none of the background functions must be computed and thus 
it is convenient to adopt ( [5.17| ). Moreover we have to notice that the zeros of background 



two-point functions and the zeros of quantum two-point functions coincide. Thus the mass 
renormalization is the same in both approaches. 

For fermions we consider the corresponding two-point function 

Tuuip) =riPLTt'ip)+ ^PR^L'ip) + PLT^u'ip) + Pn^ul'ip) , (5.18) 

where we have singled out the four different spinor structures entering in the amplitude 
^uu'{p)- An analogous equation holds for the d quark. We consider the zeros of the func- 
tion II 

■ - (r^)t(r^)-ir^) (p) , (5.19) 



Det (r^ 



where F^, and are matrices in the flavour space. The matrices in ( p.l9| ) are hermitian. 



The zeros p* of ( ^.19|) can be identifled with the masses of the quark flelds. Thus, this flxes 
the free parameters, namely the "diagonal values" of the Yukawa matrices. 

Finally, the gauge parameters ^a, C,z and have to be flxed in such a way that the 
restricted 't Hooft gauge or the Feynman gauge can be imposed. We also have to flx the 
charges cxqed^Gf and a^. This is done in the usual way |1T7[| . 

For the renormalization of the CKM matrix, there are two possible choices as was de- 
scribed at the end of Section 0. 



^^The superscript T indicates the transversal part. 
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A BFM Lagrangian and notations 

In this appendix we describe some ingredients useful for the calculation of radiative cor- 
rections in the SM within the BFM. In particular, we discuss how to derive the complete 
Lagrangian for the SM with background fields, we provide the BRST source (anti-fields) 
terms and the gauge fixing terms. 

Fields 

Here, the field content of the SM including background fields is listed. Their Faddeev-Popov 
charges, Q^n, and their UV and IR degrees^are given in the form {Qcs>n)ijv. The UV degree 
and the IR one for massless particles coincides with the mass dimension of the corresponding 
field. For massive fields the IR degree is set to two. 

• Quantum fields 

- Gauge field^3 V = {A„ Z,, G^l) with {Ol, 0?, 0?, Q\} 

- Scalar fields $ = {G^,G^,H} with {02,02,0?} 

- Ghost fields C = {c^, c^, c±, c^} with {ig, ij, ij, 

- Anti-ghost fields C = {c'^.c^ .c^.c"} with {(-1)^, (-l)i, (-1)^, (-1)^} 

- Fermion fieldsR * = {Lf , Qf , u^, d^} with {Ol , Ol , Ol , Ol , Ol } 

2 2 2 2 2 

- Nakanishi-Lautrup fields b = {6^, 6^, 6±, h"} with {0|, 0^, 0^, 0|} 

• Background fields 

- Gauge fields V = {i^, Z^, W}, 0^ with {Oj, 0?, 0?, Oj} 

- Scalar fields ^ = G±, H} with {0?, 0?, 0?} 

• External fields 

^^The UV and IR degrees are defined according to the BPHZL prescriptions given in [p^ . 
^®We use the convention — cwZf_i — sn/A^. 

-'^^The colour index is omitted. The index "i" denotes the two components of the SU{2) doublet, L^, . 

3 

In the case of the leptons the low components have the quantum numbers 0|. This is because the neutrinos 

2 

are massless. 
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BRST sources (or anti-fields) 

0* = {W*/, G^'*, c*'3, c*'±, c*'^ i/*, G^'O, LI,, L*Q^,, l\ u\ d*] with 

(-1)1 (-l)t (-2)1, (-2)t (-2)i (-l)i, (-l)i, (-l)i, (-1)1, (-1)1, 
(-1)1, (-1)1, (-1)1} 

2 2 2 

Background sources 

n ^ n^, n^, n'} with {(i)?, (i)? , (i)l, (i)?, (i)?, (i)?} 



Lagrangian and Feynman rules 

To derive the Lagrangian of the SM within the BFM we consider the invariant Lagrangian 
jC-iNviV, ^, ^, ^]{x) presented in and we replace the gauge and scalar fields by the sum 
of the quantum and the corresponding background field: Cinv[V, ^, ^, — > CiNviV + 
V, $ + $, \E', 'I'](x). It is easy to check that this procedure provides the Feynman rules of [p5|| . 
However, going beyond the one-loop level, one needs the BRST-source terms to derive the 
Feynman rules and also the WTIs and STIs. 

For the BRST-source terms we have the following equation 



BRST 



'^^ ± ( A '^w z 



=F— c c 
2 ' 



H* 



le 



2s 



w 



G+ + G+ c"- (G- + G-]c 



G° + G° 



G 



2sw 
T ie (G^ + G 



2c\ySw j 



+ G 



*,o /_ 



\2sw 



G^ ^G^\c- ^[G- ^G-\c 



2s\YC]y 



( ie 



\f2sw 

ie 



\ V2. 



LdC^ — ie 
LiiC te 



( 1 



\2s\yC\y 



^ -^w^ z 



Sw 



A ( ^ 



\2sy/c-\^ 



La 
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/ ieVud L + ■ 



i<l*u, <l*d) 



Sw 



Qdc^^ 



— g.^ 

2swCw Cw 



Sw 



2swCw Cw 



\ ^/2sw 



Z / \ 



(A.l) 



The gauge fixing terms, the Faddeev-Popov terms and the terms depending on f2 are 
easily obtained by 

I^GF + /:$n = sn{x) , (A.2) 

where is given by 

n{x) = 



J2 c"J^° 

a=A,Z,±,a 



] 



+i e ez^^^ ^ {g- G+ - G-) - e ^ztt^ — (g'' H - H G'' + 



Ciy S^ 



2cy[rS\Y 



(A.3) 



s is the BRST transformation which can be read from the Eqs. ( A.l ) and (2^). Notice 
that by taking into account the extended BRST transformations for the background fields 
(cf Eqs. ( |2.19| )) the action of s on the gauge fermion H produce new terms. An analysis of 



these new Feynman rules has been presented in |21 



Background gauge transformations 

Here we present the Background gauge transformations (BKG) for the fields and the sources. 

• BKG transformations for gauge fields 

In the next expression, Aa, A^, \± and Aa are the infinitesimal parameters of the gauge 
group SUc{S) X SUi{2) x f/y(l). Aa generates the transformations of the subgroup 
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Uq{1) of QED. 



hZ^ = d^Xz-te^{w;^X_-W-X 



Sw 



5aG^ = 5^A„-^.r'''=G;;Ae (A.4) 



BKG transformations for scalar fields 



5aG± = ±-^(i^±^G°)A±T^eG^ A^-^^^aJ , 

5aG± = + ^; ± A± T ^eG^ f A^ - 4^^A, 

ie 



5xH 



2sw 

ie 

2sw 
e 

2^ 



-G- 






-G- 






+ G- 


■A.) 


2sy/Cw 


+ G- 


-A,) 





= -^(G+A_ + G'-A+) (if + ^;)Az (A.5) 

2Si4/ ^ 

Here one must notice that the quantum scalar fields undergo a homogeneous transfor- 
mation, while the background scalar fields transform according to the inhomogeneous 
ones where the shift of the vacuum v appears. 

BKG transformations for ghost fields, for anti-ghost and Nakanishi-Lautrup 

fields 

c"^ [Xa Xz + A± c c 



Sw ^ ^ Sw 

z ■ cw 



5xc^ — ^ie 

8\c^ = — ie-^ fc^A_ + A+c~) , 
Sw ^ ' 

8xc^ = ie (c+A_ + A+c~) , 



5ac" = -gsr'^c'X, (A.6) 
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The BKG transformations for the Nakanishi-Lautrup fields and for the anti-ghost fields 
are identical with the previous ones after substituting the ghost fields with the other 
fields. 



• BKG transformations for the fermion fields 



Sxf" 



le 



\/2sw 

ie 



LdX+ - ie QuXa - 

L„A_ - ie Qd><A + 

-ieQi (Xa + -^Xz) , 
ieVud 



ieV 



Sw 



Cw 

-qdX+ - ie 



ud 



?„A_ - ie 



QuXa — 

QdXA + 



ASwCw Cw 



1 n ^^\\' 

■z — I '^Z 

ZSwCw Cw / 

ASwCw Cw J 



Ld , 



lu + i^Qs {r^qu) Xa , 
Qd + ^9s {r^qd) A„ , 



V^Sw 

-ieQf (a^ + ^Az) /« + '-g, (r"/^) A„ (A.7) 

Vud is the CKM matrix, f = u,d and L = {u, e). 

• BKG transformations for the BRST sources 

The BKG transformations of BRST sources correspond to the gauge transformations 
of the corresponding quantum gauge field according to their specific representations. 



Linearized Slavnov- Taylor operator and functional Taylor operator 

The linearized Slavnov-Taylor operator for a generic functional T is given by 



where 



St{^) = J d'^x{swdfj,cz + cwdf^CA) (^sw^ + cw^^^ 



+ E h^— + {T,T) + {:F,T), 



a=A,Z,±,a 



(A.8) 



SX 5Y 6X SY 5X 5Y 6X SY 5X SY 



5X 5Y 5X 5Y 



6X 6Y 6X 6Y 6X 6Y 
+ s^s^ + 5^5^ ^ 5G*^5G^ ^ SG*^SG^ ^ 5H*5H 
I SX SY 



+ h.c. 



(A.9) 



Since S{T) = 0, the operator St is nilpotent. We also introduce the tree-level linearized 
operator Sq which is equivalent to »Sro where Fq is the tree level action. 
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The Taylor operator on the functional T is defined as follows. One first considers 
the relevant amplitude which results from functional derivatives w.r.t. fields denoted by 
subscripts ^ (i}i{pi)(f)2{p2)...(i>m{pm) with YlJLiPj = 0. Then the Taylor expansion in the inde- 
pendent momenta up to degree 6 acts formally as 

oo „ m 

T'r = E I[<^'p^Mp^)s\EUP^K,■..,p,. ^MPMP.y-MpM^ ^ , ■ (A.io) 



m=l i=l 



A remarkable property of is that T^^T^'^ = with S = min{5i,52}- Note that the 
Taylor operator is scale invariant, however, it does not commute with spontaneous symmetry 
breaking. If massless fields are present IR problems can occur. 



B An example for the derivation of STIs 

In this appendix we explicitly derive the STIs for the amplitude involving a scalar matter 
field $ and two gauge fields V^* and Thereby we follow the general rules derived in 
Section pTTI. The amplitude for this process is computed in terms of the following irreducible 



Green functions: the two-point functions Fvmt/^, Ft/m^ (i,j = 1,2) and which repre- 
sent corrections to external legs and the three-point function FyMy^i.^. The renormalization 
program also requires the correct definitions of other Green functions which are absent in 
the physical amplitude which, however, arise in the analysis of the STIs. 

Let us in a first step discuss the two point functions Tyt^yv. The ghost field, associated 
with the gauge boson Vf, is denoted by Cj. According to rule ^ we have to differentiate 



Eq. w.r.t. Ci and V"/: 



5ci{-p)5V^{p) 



<t>=o 



+ ^'^^yip (p)^v£v- (p) + {pW^v- ij>) 

= 0, (B.l) 

where in the second line a summation over k is understood. Notice that for simplicity in the 
second line the gauge eigenstates notation is used. The subscript "0 = 0" means that after 
differentiation, all fields are set to zero. $ and V represent all the scalar and gauge fields, 
respectively. H denotes the Higgs field and $*, V* and H* are the corresponding anti-fields. 
In the following, the summation over the scalar fields is omitted. If one does not use normal 
ordering in the perturbative expansion, one gets an additional tadpole term on the r.h.s. of 
( [B.l|) , namely ]rc^H*v{^jP)^H{^)- However, this contribution is easily removed by using the 
renormalization condition F/^ (0) = which implements the spontaneous symmetry breaking 
at the quantum level. The linear terms in Eq. (|B.1D take into account the contributions 



arising from the abelian gauge field. The quantity 6ij represents the Kronecker delta carrying 
gauge field indices. Notice that no charged gauge boson is involved in these terms. 
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As mentioned below Eq. ( |2.6| ), formula ( |B.1| ) is valid to all orders. We want to stress 



that at tree level, JTyly^ can be read off the invariant Lagrangian ||T^, |T^, |T8[. 



Let us now continue the construction of the closed set of STIs for the example under 
consideration. Eq. ( |B.1| ) contains new Green functions, namely the two-point function for 
the mixing between the gauge and the scalar fields, r$yi'(p), and the ones involving one 
ghost and one anti- field, namely Tci^*{p) and ^av* ip)- According to rule |1], these Green 
functions are not vanishing and they are related to the two-point functions of ghost and 
anti-ghost through the ghost equations which reads in the 't Hooft gauge fixing ||^, [TB|, PT| : 

^c,-cAp) = (p) + OM,,<i.r,^$.(p) . (B.2) 



Mj is defined below (2.25). This equation is independent from the STIs and it implements 



the equation of motion for the ghost fields. It is a consequence of the choice of the gauge 
fixing. 

In order to get an identity involving Wq>yi'{p) we can derive a new independent STI by 
differentiating (|2.6|) w.r.t. Cj and 



5cj{-p)5^p) 



<j>=0 



0. 



(B.3) 



In this equation, only one new Green function emerges, namely r$/$(p). As it does not 
contain any gauge or ghost fields the process of finding a closed set of STIs is completed. 
There is no independent STI which can give us new information on the Green functions. 
Notice that also the Green functions involving ghost fields, Tcj^*{p) and ^cjV* {p), are the 
same as before. 

Let us now consider the three-point function Ty^y^^^. Replacing one of the gauge bosons 
by its corresponding ghost field leads to (omitting the summation over the Goldstone fields) 



(535(r) 



<f>=0 



SCi{p3)SV^ {Pl)6<l>{p2 

- iP3 {swSiZ + CwSiA) {sw'S^^^yj^{pi,p2) + CwT^^y,^{pi,p2 
+ ^'=^V';p(-P3)ir^^^^.^(pi,P2) + T^^y,y,{p2,Pl)Ty^^^{p2) 

+ rc,v;*^*(pi,P2)irv'^,,y,-(pi) + r,^$'*(-P3)ir^/yj$(pi,P2) 

+ r^,$'.yj(P2,Pi)ir$'$(P2) + T^^^u^{pi,P2)T^,yj{pi) 
= 0. 



(B.4) 



Besides the two-point functions already discussed above, new Green functions arise: T^,yj^, 
ffj'.vj ciiid The first one satisfies a new independent STI because of the presence 
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of only one gauge field. In fact, differentiation w.r.t. and $ leads to the identity 



0=0 



- ips {swSiZ + CwSiA) {sw'^Zt^'S>'-S>{Pl,P2) + cvyr^M<i,'$(pi,P2)) 



k,p 



+ T^^v*^^{pi,P2Wv^,,'s>'{pi) + ir,^$"*(-P3)ir$//$/$(pi,p2) 

0. 



(B.5) 



The summation over the Goldstone fields is understood. In this identity, only the new 
Green function r$//$/$ occurs which is not constrained by an independent STL Concerning 
the Green functions with only external quantum fields the system is now complete. However, 



we are left with the ones involving also ghost fields. In the case of F 



w.r.t. Ci, Cj and ^4*^ and get 



we differentiate 



<t>=o 



6ci{ps)6cj{pi)6VC^p{p2 

- iPs {sw^iz + cwSiA) {sw'S^zi^c,vi^{Pi,P2) + cw'ir 

+ ip'l {swSjZ + CwSja) [sw'SI'zt^c,V*^{p3,P2) + Cw'S:'At^c,V*^{p3,P2 



A^c,V,*JPl,P2 



+ ^c,v:i-Pi)T 



+ re^$,(-pi)r$c,v;*/P3,P2) + rc,c,cr(pi'P2)irc,v;*/P2) 



0. 



{Pl,P2) 



(B.6) 



An analogous equation is obtained for T^^^uyj- In ( p.6|) c* is the anti-field of Cj. At this 
point a short comment is in order. As is well known, the BRST transformations are non- 
linear and local. At the quantum level, they receive radiative corrections. Moreover, if 
the regularization breaks the symmetries, the correct non-linear transformation rules at the 
quantum level can be obtained only by adjusting the finite counterterms which are needed 
for the STI (p.6| ). Note that they — in contrast to (|B.4|) and (|B.5|) — contain derivatives 
w.r.t. anti-fields. 

In ( p.6| ) the new function Wc-c^c* which does not contain gauge fields occurs. According 
to rule HI we get the following identity: 



0=0 



SCi{p4)6Cj{pi)6Ck ip2 ) 5 C*^ {P3 ] 

- ipt {sw^iZ + Cw^iA) (sw^ Z^^c,CkC*SPl^P^2^P3) + Cw^r^Mc.XfcC^ (Pl, P2, Ps^ 
+ ipt {swSjZ + CwSja) {sw'S^zi^c,CkC*^{P4,P2,P3) + Cvp-r^Mc.cfcC^ (P4,P2,P3; 
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- ipi^ {swh,Z + Cwh,A) {sw^ Zf^c,c,c*^{P4,Pl,P3) + Cw'S^Ai-c,c,c*^{P4,Pl,P3)) 

+ ^c^V;S-P^)^V^'^c,cjc*^{P4,Pl,P3) + irc,<I.*(-p2)ir$c.c,cJ„(P4,Pl,P3) 

+ ^c,V*^i-Pl)'S^V^^^^CkC,c-t^iP2,P4,P3) + rc^**(-Pl)r$CfcC,cJ„(p2,P4,P3) 

+ rc,y,*^(-p4)ryi,Mc,CfcC^(Pl,P2,P3) + ^ c,'S>* {-Pi)'^ 'S>c,CkC*^{pi, P2, Ps) 

+ ^c,c^c;{PuP2+P3)'S^cic,c*^{P2,P3) 

+ ^c,c,c;iP3,P2+PlWctc,c'JP4,P3) 

+ ^c,c,c;{p4,Pl +P3)Tc,c,c'^{pi,P3) 

= 0, (B.7) 

which can be considered as the extension of the Jacobi identity to the quantum level. 

From rule g it also follows that derivatives w.r.t. CjCj V^^VJ"', CiCj$$*, CjCjV^^$* have to 
be considered. The corresponding identities do not involve any new Green function which 
requires a subtraction, however, they provide new constraints^. For brevity we will not list 
them explicitly. 

Now all symmetry constraints for the amplitude and the related Green functions have 
been derived. In a second step, the same procedure has to be applied to the sub-divergences. 
Again the independent irreducible contributions have to be extracted and a closed set of 
independent STIs has to be derived. At that point, one is able to study the counterterms 
and the breaking terms involved in the set of identities. 



C Triangular organization of counterterms 

As is clear from the applications discussed in Sections § and ^ and from the theoretical anal- 
ysis of Section |3.2| , the triangular organization of functional identities and of counterterms 
is very important. With triangular structure of functional identities we mean the possibility 
to organize the set of functional identities into a hierarchical structure in such a way that 
we can restore the identities one after the other without spoiling those which are already 
recovered. As an example, we recall the 6 — >■ 57 calculation. We have seen that we can 
organize the counterterms in such a way that in a first step the WTIs and in a second one 
the STIs can be restored. This can only be achieved if the fixing of the STIs does not de- 
stroy the already restored WTIs. This means that the counterterms needed to restore the 
STIs must be invariant under the action of the WTIs, i.e. they must be background gauge 
invariant. Clearly this is only possible if the breaking terms to the STIs are background 
gauge invariant. In Section |3.2| , by using the consistency conditions and repairing the WTI 

^''Here we would like to underline that the problem of the complete set of identities involved is a pure 
algebraic problem. This is equivalent to derive the well-known descent equations in the space of local 
functionals If anti- fields are required, the relevant constraints can easily be derived from the descent 
equation formulated in the Batalin-Vilkovisky anti- field formalism . As an example, consider the ampli- 
tude for VV —> cV. According to rule 1 the STI for cV produces cV* . The latter requires the identity for 
cV* ^ c'^ V* which finally produces c^c*. 
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by suitable non-invariant counterterms, we have shown that the new breaking terms 

A'5 = A', - 5^ (/ dS/:r/,) (c.i) 

are indeed background gauge invariant W"(a) (A's) = 0. 

We mention here that this fact can be extended to the complete set of functional identities. 
We recall that the SM in the BFM ||2l| is completely defined in terms of the following 



functional identities (up to free arbitrary constant parameters) 



The Nakanishi-Lautrup identities ( p.7[) , which implement the gauge fixing conditions 
to all orders, 



the Faddeev-Popov equations (see l2l|, 
the Abelian Anti-ghost Equation (in the case of BFM ||21||), 
the WTI given by Eq. ( p. 241) for the background gauge invariance and 
the STI given by Eq. ( |2.6| ) for the ERST symmetry. 



We specialize now the argument to each subspace of counterterms and we show that it is 
indeed possible to restore the WTIs, respectively, the STIs by starting from the identities for 
two-point functions, then for three-point functions and, finally, to fix the counterterms for 
four-point functions. Clearly, this is a considerable simplification for practical calculations. 
It helps us to restrict the set of identities involved into a specific calculation and, as a 
consequence, it restricts the number of counterterms which must be effectively restored. 

We first focus on the WTI. The main problem in dealing with the WTI is that the operator 
W(\-) is not homogeneous in the fields. Since this operator implements the local background 
gauge symmetry, the transformations of the gauge fields contain an inhomogeneous term 
which does not depend on the gauge fields itself. Furthermore, the transformations of the 
scalar fields contain the shift of the Higgs fields in order to take into account the spontaneous 



symmetry breaking mechanism. Following the extended anti-field formalism []T9| we promote 
the infinitesimal parameter A to a local dimensionless Grassman parameters uj in the adjoint 
representation of the gauge group which transforms asuj^uj + ujf\uj and we decompose 

+ (a; A ^),^ = >VS)+>V('i), (C.2) 
where W^'^) contains the transformations of uj. This definition implies that 

(Wj))' = 0, (W('i))' = 0, {wS),W('i)} = 0. (C.3) 

On the other side, the action F and the breaking term A'^y(ti;) can be decomposed into am- 
plitudes with n external legs0 F = '}2n>i ^^ct respectively A'^^^cu) = Y.n=i ^nwi^) where 

^^This means that r„ is generically an n-point function and are coefficients of monomials with n fields 
(comprehensive of uj). 
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is the upper bound. In the case of renormahzable theories we have = 4. Therefore, the 
contribution with the highest content of fields, namely n = — 1, A^, to Eqs. (|3.11|) are 

w'i(rjv) = A^^i,^ - >v'° (r^) 
(r^v+i) = A'j,^w - (r^+i) , (C.4) 

where W"°'^ are the matrices acting on each single monomials of r„. The breaking terms 
satisfy 

>v'i(a^,^) = 

>V'° (a;.,^) + W'^ (a^_i,^) = 0. (C.5) 

Notice that >V'° (r^r+i) and W'^ {S^n) are combinations of amplitudes with A^ external fields 
whose finite parts can be adjusted by counterterms with monomials with A^ fields. On the 
other side, W"^ (fAf) are combinations of amplitudes with A^ — 1 external fields whose finite 
parts are already fixed by the WTI for breaking terms with 0, . . . , A^ — 2 external fields. 
Finally, W'^ (r^v+i) is expressed in terms of A^ + 1 external fields amplitudes which cannot 
be modified by adjusting overall local counterterms. 

Therefore the two systems of Eqs. ( |(J.4| ) must be solved simultaneously and this can be 
done only if the system is redundant. To show that we apply the operator W"^ from the 
right on the first equation and the operator on the second and by using the consistency 
conditions ( |C5[ ), we immediately find that a non-trivial combination of the two system of 
equations exists which implies the redundancy of the two system. The consistency conditions 
help us to disentangle the independent equations. 

For the STI the problem is a little more involved. The functional operator iSg is not 
linear, and decomposing it w.r.t. the power of the fields we get a mixing between terms 
coming from the functional operator and those coming from the action F. 

We introduce the notation iSq^ (A^ = 1,2,3,4) to denote the contribution linear, 
quadratic, cubic and quartic to the non-linear Slavnov- Taylor operator and, as above, we 
indicate with A'j^ grpj {N = 1,2,3,4,5) the breaking terms of the STI. By expanding the 
equation Sq{T) = A'^j^j and the consistency conditions S^^Agrpj) = in powers of fields, we 
have the following equations 

'^0,2^2 = A2,5T/ 5 
'^0,2^2, 5T/ = ) 

'^0,3^2 + Sq 2^3 = A'^^gj^i , 
'^0,3^2, 5T/ + '5o,2^3,ST/ = 5 

'5o^4F2 + ^o^gFs + ^0,3^4 = A'^,STI 5 
^0,4^2,3X1 + '^0,3^3,ST7 + '^0,2^4,ST/ = ) 

«5o,4F3 + 5o 3F4 + 3^ 2^^ = A'^^gj.j , 
A' _i_ A' _i_ A' — n 

Rearranging the last two equations, we can immediately see that they must compatible 
by means of the consistency conditions otherwise they cannot be solved in terms of the 



59 



remaining free parameter Notice that the Green functions Fs are superficially convergent 
and therefore they cannot be fixed by the previous equations. We apply Sq from the left on 
both equations and use the commutation properties 

(5^J' = o, (5^,3)' = 0, (5^,4)' = 0, 

in combination with the consistency conditions. Then one checks the compatibility between 
the two equations. This completes the proof of the triangular structure of the counterterms 
for the STIs. 

References 

[1] C.N. Yang and R.L. Mills, Phys. Rev. D 96 (1954) 191; T.D.Lee and C.N. Yang, 
Phys. Rev. D 98 (1955) 101; S.Weinberg, Phys. Rev. D 19 (1967) 1264; A. Salam, 
in Elementary Particle Physics, ed. N. Svartholm (Almqvist and Wiksells, Stockoholm, 
1968); S.L. Glashow, Nucl. Phys. B 22 (1961) 579; A.Salam and J. Ward, Phys. Lett. B 
13 (1964) 168; J. Schwinger, Ann. Phys. (NY) 2 (1957) 407. 

[2] CO. Bolhni and J.J. Giambiagi, Phys. Lett. B 40 (1972) 566; G. 't Hooft Nucl. Phys. B 
61 (1973) 455; G. 't Hooft, Nucl. Phys. B 62 (1973) 444; D.A. Akyeampong and 
R. Delbourgo, Nuovo Cim. 17A (1973) 578; J.C. Collins, Nucl. Phys. B 80 (1974) 341; 
W.J. Marciano, Nucl. Phys. B 84 (1975) 132; Nucl. Phys. B 88 (1975) 86; T. Marinucci 
and M. Tonin, Nuovo Cim. 31A (1976) 381. 

[3] G. 't Hooft and M. Veltman, Nucl. Phys. B 44 (1972) 189. 

[4] P. Breitenlohner and D. Maison, Comm. Math. Phys. 52 (1977) 11; Comm. Math. 
Phys. 52 (1977) 39; Comm. Math. Phys. 52 (1977) 55. 

[5] M. Ciuchini, E. Franco, G. Martinelli, L. Reina, and L. Silvestrini, Phys. Lett. B 316 

(1993) 127; M. Ciuchini, E. Franco, L. Reina, and L. Silvestrini, Nucl. Phys. B 421 

(1994) 41. 

[6] R. Harlander and M. Steinhauser, Prog. Part. Nucl. Phys. 43 (1999) 167. 

[7] H. Epstein and V. Glascr, Ann. Inst. Poincare 29 (1973) 211, Adiabatic Limit in Per- 
turbation Theory, in G. Vclo, A.S. Wightman (cds.), Renormalization Theory; D. Rcidcl 
Publishing Company, Dordrecht 1976, 193; R. Stora, Differential Algebras, ETH-lectures 
(1993) unpublished; T. Hurth, Ann. Phys. (NY) 244 (1995) 340; T. Hurth and K. Sk- 
enderis, Nucl. Phys. B 541 (1999) 566. 

[8] Y.M.P. Lam, Phys. Rev. D 6 (1972) 2145; Phys. Rev. D 7 (1973) 2943; W. Zimmer- 
mann. Comm. Math. Phys. 39 (1974) 81; Local Operator Products and Renormaliza- 
tion in Quantum Field Theory in 1970 Brandeis University Summer Institute Lectures, 



60 



Cambridge, Mass. M.I.T. Press and references therein; J. Lowenstein, Comm. Math. 
Phys. 24 (1971) 1; J.H. Lowenstein and B. Schroer, Phys. Rev. D 7 (1975) 1929. 

[9] C. Becchi, A. Rouet, and R. Stora, Comm. Math. Phys. 42 (1975) 127; Ann. Phys. (NY) 
98 (1976) 287; I.V. Tyutin, Lebenev Institute preprint N39 (1975); see also C.Becchi, 
Lectures on Renormalization of Gauge Theories in Relativity, group and topology II, 
Les Houches 1983. 

[10] O. Piguet and A. Rouet, Phys. Reports 76 (1981) 1 and references therein; O. Piguet 
and S.P. Sorella, Algebraic Renormalization Lecture Notes in Physics Monographs, p 
28, Springer- Verlag Berlin Heidelberg, 1995 and references therein. 

[11] G. 't Hooft, Nucl. Phys. B 33 (1971) 436; H. Kluberg-Stern and J. Zuber, Phys. Rev. D 
12 (1975) 467; Phys. Rev. D 12 (1975) 482; L.F. Abbott, Nucl. Phys. B 185 (1981) 
189; S. Ichimose and M. Omote, Nucl. Phys. B 203 (1982) 221; D.M. Capper and 
A. MacLean, Nucl. Phys. B 203 (1982) 413; D.G. Boulware, Phys. Rev. D 12 (1981) 
389. 

[12] W.A. Bardeen, A.J. Buras, D.W. Duke, and T. Muta, Phys. Rev. D 18 (1978) 3998. 

[13] J.H. Lowenstein and W. Zimmermann, Comm. Math. Phys. 44 (1975) 73; Comm. Math. 
Phys. 46 (1976) 105; Nucl. Phys. B 86 (1975) 77; J.H. Lowenstein, Comm. Math. 
Phys. 47 (1976) 53; Nucl. Phys. B 96 (1975) 189; T.E. Clark and J.H. Lowenstein, 
Nucl. Phys. B 113 (1975) 109. 

[14] J. Zinn- Justin, in Trends in Elementary Particle Theory, Lecture notes in Physics vol. 
37, Springer Verlag, Berlin, 1975, eds. H. RoUnik and K. Dietz; B. Lee, in ''^Methods 
in Field Theory, Proceedings Les Houches session 28", 1975, eds. R. Balian and J. 
Zinn- Justin, North-Holland, Amsterdam, 1976. 

[15] W. Zimmermann, Local Operator Products and Renormalization in Quantum Field The- 
ory in 1970 Brandeis University Summer Institute Lectures, Cambridge, Mass. M.I.T. 
Press. 

[16] G.Bandelloni, C.Becchi, A.Blasi, and R.CoUina, Ann. Inst. Henry PoincareXXMlll, 3 
(1978) 225, 285; C.Becchi, A. Rouet, and R.Stora, Renormalizable Theories with Sym- 
metry Breaking in Field Theory, Quantization and Statistical Physics, Ed. E. Tirapegui; 
G. Bandelloni, C. Becchi, A. Blasi, and R. Colhna, Comm. Math. Phys. 71 (1980) 239. 

[17] M. Bohm, H. Spiesberger, and W. HoUik, Fortschr. Phys. 34 (1986) 687. 

[18] E. Kraus, Ann. Phys. (NY) 262 (1998) 155. 

[19] 1. A. Batahn and G. A. Vilkovisky, Phys. Lett. B 102 (1981) 27, Phys. Rev. D 28 
(1983) 2567; M. Henneaux and C. Teitelboim, Quantization of Cauge Systems, Prince- 
ton University Press, Princeton, 1992; J. Gomis, J. Paris, and S. Samuel, Phys. Reports 
259 (1995) 1; G.Barnich, F. Brandt, and M. Henneaux, Comm. Math. Phys. 174 (1995) 



61 



93; Comm. Math. Phys. 174 (1995) 57; Phys. Lett. B 346 (1995) 81; G. Barnich and 
M. Henneaux Phys. Rev. Lett. 72 (1994) 1055; F. Brandt, M. Henneaux, and A. Wilch, 
Nucl. Phys. B 510 (1998) 640. 

[20] K. Aoki, Z. Hioki, R. Kawabe, M. Konuma, and T. Muta, Prog. Theor. Phys. Suppl. 
73 1982 1. 

[21] P.A. Grassi, Nucl. Phys. B 462 (1996) 524; Nucl. Phys. B 537 (1999) 527; Report No.: 
MPI/PhT/98-57, submitted to Nucl. Phys. B. 

[22] C. Itzykson and J.C. Zuber, Quantum Field Theory, McGraw-Hill, 1980. 

[23] E.R. Spccr, "Dimensional And Analytic Renormalization," In *Erice 1975, Proceedings, 
Renormalization Theory*, Dordrecht 1976, 25-93 and references therein. 

[24] P. Breitenlohner and H. Mitter, Nuovo Cim. lOA (1971) 655; Nucl. Phys. B 7 (1968) 
443. 

[25] A. Denner, G. Weiglein, and S. Dittmaier, Phys. Lett. B 333 (1994) 420; Nucl. Phys. B 
440 (1995) 95. 

[26] J. Papavassiliou, Phys. Rev. D 41 (1990) 3179; J. Papavassiliou and A. Pilaftsis, 
Phys. Rev. D 54 (1996) 5315. 

[27] R. Ferrari, A. Le Yaouanc, L. Oliver, and J.C. Raynal, Phys. Rev. D 52 (1995) 3036. 

[28] J.C. Romao and A. Barroso, Nucl. Phys. B 272 (1986) 693; J.M. Soares and A. Barroso, 
Phys. Rev. D 39 (1973) 1989. 

[29] P. Gambino, P.A. Grassi, and F. Madricardo, Phys. Lett. B 454 (99) 98. 

[30] W. Marciano and A. Sirlin, Nucl. Phys. B 93 (1975) 303; A. Denner and T. Sack, 
Nucl. Phys. B 347 (1990) 203; B. Kniehl and A. Pilaftsis, Nucl. Phys. B 474 (1996) 
286. 

[31] A. Strumia, Nucl. Phys. B 532 (1998) 28. 
[32] S. Weinberg, Phys. Rev. D 118 (1960) 838. 

[33] W.A. Bardeen, Phys. Rev. 184 (1969) 1848; J. Dixon, Cohomology and Renormalization 
of Gauge Theories (1976,1977) unpublished. Comm. Math. Phys. 139 (1991) 

495; B. Zumino, Chiral Anomalies and Differential Geometry, in Relativity, Groups 
and Topology II, ed. B.S. DeWitt and R. Stora, North Holland, Amsterdam (1984); R. 
Stora, Algebraic Structure and Topological Origin of Anomalies, in Progress in Gauge 
Field Theory, ed. 't Hooft ct al, Plenum Press, New York (1984); F. Brandt, N. Dragon, 
and M. Kreuzer, Phys. Lett. B 231 (1989) 263; Nucl. Phys. B 332 (1990) 224, 250; M. 
Dubois- Violette, M. Henneaux, M. Talon, and CM. Viallet, Phys. Lett. B 289 (1992) 
361. 



62 



[34] R. Ferrari and P.A. Grassi, |hep-tJi/9807191| , to appear in Phys. Rev. D; R. Ferrari, 
P.A. Grassi and A. Quadri, |hep-th/ 990519^ . 

[35] M. Veltman, Nucl. Phys. B 123 (1977) 89. 

[36] M.A. Shifman, A.I. Vainshtein, M.B. Voloshin, and V.I. Zakharov, Sov. J. Nucl. 
Phys. 30 (1979) 711. 

[37] A. Djouadi, M. Spira, and RM. Zerwas, Phys. Lett. B 264 (1991) 440. 

[38] C. Greub and T. Hurth, Nucl. Phys. Suppl. B 74 (1999) 247 and references therein. 

[39] A. Czarnecki and W.J. Marciano, Phys. Rev. Lett. 81 (1998) 277. 

[40] J.F. Donoghue, Phys. Rev. D 19 (1979) 2772. 



63 



